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A  "regenerative  load  system"  is  proposed,  replacing  the  conventional  load  for 
the  long-term  testing  of  power  equipment  and  motor-generator  set.  Basically,  this 
system  consists  of  two  parts:  (1)  impedance  emulator  and  (2)  inverter  with  power 
filtering  function. 

The  proposed  system  has  three  significant  merits:  economic,  safety  and  power 
quality.  The  first  merit  involves  the  recycling  of  the  electrical  energy,  resulting  in  the 
energy  consumption  less  than  20%  of  the  conventional  load.  If  the  soft-switching 
technique  is  employed  for  the  converter,  the  total  energy  dissipation  can  be  even  lower 
than  10%.  Second,  safety  is  a  consideration.  Due  to  the  low  energy  dissipation,  no 
appreciable  Joule  heating  is  observed.  The  third,  power  quality,  is  equally  important. 
Since  the  current  from  the  ac  line  is  shaped  into  a  pure  sinusoidal  waveform  with 


power  factor  equal  to  one,  problems  such  as  harmonic  currents  and  EMI  will  be  greatly 
reduced.  Hence  the  power  quality  of  the  utility  source  can  be  improved. 

In  order  to  accomplish  the  analysis  and  design  of  the  system,  the  following 
contributions  have  been  made.  For  the  active  power  filter  design,  a  new  control 
method,  using  an  integrator  and  sampling  technique,  is  proposed  to  simplify  the 
calculation  algorithm  for  the  real  fundamental  component  of  load  current.  In  addition, 
a  new  simple  control  scheme,  based  on  the  energy  balance  concept,  is  proposed  to 
control  the  voltage  of  an  energy-storage  capacitor.  The  advantages  of  this  active 
power  filter  (APF)  circuit  are  simplicity  of  control  circuits,  low  cost  (a  smaller  energy- 
storage  capacitor)  and  good  transient  response  (less  than  two  cycles  of  60  Hz).  The 
feasibility  of  this  theory  is  verified  by  using  PSpice  simulation  and  experimental 
results. 

Secondly,  in  the  design  of  the  system,  a  memory  circuit  is  used  to  calculate  and 
memorize  the  utility  source  current.  This  design  provides  an  alternative  method  to  the 
design  of  APF  control  circuits.  Third,  we  have  found  that  the  utility  source  inductance 
may  cause  the  APF  system  to  become  unstable.  After  the  addition  of  the  damping 
circuit,  the  system  can  be  stable.  However,  even  with  the  addition  of  the  damping 
circuit,  the  APF  or  regenerative  load  system,  will  still  be  unstable  if  the  gain  of  the 
current  controller  is  greater  than  one.  Therefore,  for  the  stability  of  the  system  with 
utility  source  inductance,  the  following  conditions  must  be  met:  (1)  addition  of  a 
damping  circuit  and  (2)  the  current  controller  gain  must  be  kept  smaller  than  one.  This 


vi 


phenomenon  has  been  verified  by  Matlab  control  analysis  and  PSpice  simulation  and 
has  never  been  noted  in  the  literature. 

Fourth,  employment  of  the  equivalent  impedance  concept  separates  the 
regenerative  load  system  into  two  equivalent  subsystems  and  establishes  the  "input 
impedance  specification"  of  the  DUT  (Device  Under  Test)  for  the  stability  of  the 
"regenerative  load  system."  Fifth,  for  the  improvement  of  stability  margin  of  the 
system,  a  new  "two-damping-circuits"  concept  is  proposed.  The  addition  of  a  second 
high  frequency  damping  circuit  can  further  improve  the  stability  margin  of  system. 
This  "two-damping-circuits"  concept  can  be  applied  also  to  other  circuit  design. 

Sixth,  and  lastly,  when  using  the  active  damping  circuit  to  replace  the 
conventional  passive  damping  circuit,  the  overall  system  efficiency  is  improved.  Also, 
a  new  "reactive-current-compensation"  concept  is  developed  to  compensate  the 
reactive  component  absorbed  by  the  damping  circuits.  Hence  the  power  factor  of  the 
utility  source  current  can  be  increased  to  nearly  one,  even  in  the  presence  of  the 
damping  circuits. 
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CHAPTER  1 
INTRODUCTION 

If  we  test  power  device  or  equipment  where  a  conventional  load  is  used,  a  large 
amount  of  energy  will  be  dissipated  and  Joule  heat  produced.  This  dissertation  develops 
a  "regenerative  load  system"  for  long  term  testing  of  "power  equipment"  and  "motor- 
generator  set"  to  achieve  energy  conservation.  At  the  same  time  this  system  can  improve 
the  power  quality  of  the  utility  source.  Basically  it  consists  of  two  parts:  (1)  an 
impedance  emulator  and  (2)  an  inverter  with  power  filtering  function,  as  shown  in  Fig. 
1.1.  Two  topics  in  power  electronics  are  closely  related  to  this  study:  "power  factor 
correction"  [1-23]  and  "active  power  filter"  [24-56].  The  following  discusses  some  of  the 
characteristics  and  applications. 

1 . 1  Impedance  Emulator  ( related  to  power  factor  correction ) 
For  applications  in  which  line  pollution  in  terms  of  harmonic  content  and 
displacement  power  factor  of  the  input  current  is  of  main  concern,  ac/dc  converters 
featuring  almost  unity  power  factor  are  required.  These  high  quality  rectifiers,  also  called 
"power  factor  pre-regulators"  (PFPs)  replace  the  conventional  capacitive-filter  rectifiers, 
which  have  the  disadvantages  of  absorbing  high  peak  currents  from  the  utility  line.  The 
goal   of  a  high   quality   rectifier   is   to   emulate   a   resistive   load,  achieving 
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Fig.  1  1  The  regenerative  load  system  and  the  internal  structural  block  diagram. 

(a)  The  regenerative  load  system  and  (b)  The  internal  structural  block  diagram  of  the 

regenerative  load  system. 


theoretically  the  unity  power  factor,  even  in  the  presence  of  distorted  line  voltage.  In 
order  to  do  that,  PFPs  (power  factor  pre-regulators)  must  be  able  to  shape  the  input 
current  in  such  a  way  that  it  represents  a  scaled  replica  of  the  line  voltage.  Among  the 
active  methods  for  input  current  shaping  based  on  switching  dc/dc  converters,  those 
employing  boost  topology  are  mostly  applied  [1-3].  Moreover,  this  resistor  emulator 
concept  has  also  been  extended  to  the  "variable  active-passive  reactance"  (VAPAR), 
which  was  originally  proposed  in  references  [24]  to  [26],  and  is  one  of  the  new 
applications  of  the  bilateral  converter.  The  VAPAR  works  as  a  variable  passive  reactive 
element  (or  a  kind  of  virtual  reactance)  including  the  transient  state.  Realization 
conditions  and  experimental  results  were  introduced  by  Funato  and  Kawamura  [24].  In 
these  analyses  and  experiments,  however,  a  dc  voltage  source  was  used  as  an  energy 
storage  element.  If  a  dc  voltage  source  is  replaced  with  a  capacitor,  the  VAPAR  can  be 
considered  as  a  complete  two-terminal  element.  In  Chapter  3,  the  VAPAR  concept  is 
employed  to  implement  the  "impedance  emulator." 

1.2  Inverter  with  Active  Power  Filtering  Function  (related  to  active  power  filter) 
Power  electronics  are  widely  used  in  industrial  equipment,  such  as  frequency 
changer,  motor-drive  system,  etc.  Such  equipment  presents  a  nonlinear  impedance  to  the 
utility,  generating  large  harmonic  currents  with  adverse  effects,  such  as  low  power  factor, 
low  efficiency,  and  destruction  of  other  equipment  (e.g.  the  power  capacitor  can  be 
damaged  by  the  resonant  overvoltage,  etc.).  Also,  EMI  interferes  with  instruments  and 
communication  equipment.  Therefore,  the  active  power  filter  was  used  to  reduce  the 
above  pollution  problems  and  to  compensate  the  reactive  power  of  the  load 
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simultaneously.  Since  the  active  power  filter  must  provide  the  harmonic  and  reactive 
power,  it  is  basically  a  bilateral  converter.  Chapter  2  proposes  the  active  power  filter  with 
"reduced  energy-storage  capacitor."  Subsequently,  in  Chapter  3,  this  APF  is  re-modified 
as  an  "inverter  with  power  filtering  function"  to  become  a  part  of  the  "regenerative  load 
system." 

1 .3  Regenerative  Load  System 
In  Chapter  3,  the  above  two  parts  are  combined  together  to  form  a  regenerative 
load  system,  used  for  long  term  testing  of  the  power  equipment  and  motor-generator  set. 

The  proposed  system,  shown  in  Fig.  1.1(a),  consists  of  two  parts:  (1)  an 
impedance  emulator  and  (2)  an  inverter  with  power  filtering  function.  The  internal 
structural  block  diagram  is  shown  in  Fig.  1.1(b).  In  the  first,  the  impedance  emulator 

V, 

draws  currents  from  the  DUT,  according  to  the  equation  I L  =  — — .  This  current  related 

power  can  be  used  to  charge  the  "energy-storage  capacitor."  Therefore,  the  input  of  this 
converter  can  be  viewed  as  an  equivalent  load  impedance  of  value  "ZL".  The  energy 
stored  in  the  capacitor  is  inverted  into  ac  current  and  is  fed  back  into  the  input  of  the  DUT 
by  the  second  part  of  the  system,  i.e.,  "inverter  with  active  power  filtering  function."  In 
addition,  this  inverter  can  also  provide  the  compensation  current  to  shape  the  ac  line 
current  is  into  a  pure  sinusoidal  waveform  with  power  factor  P.F.=1.0. 

The  proposed  system  has  many  merits.  First,  this  system  is  economical  because 
of  the  recycling  of  the  electrical  energy.  The  energy  consumption  will  be  less  than  20% 
of  the  conventional  load.  If  the  "soft-switching"  scheme  is  employed  for  the  converters. 


the  energy  dissipation  can  be  even  lower  than  10%.  Secondly,  safety  is  a  consideration. 
Due  to  the  low  energy  dissipation,  no  appreciable  Joule  heating  is  observed.  Lastly,  since 
the  current  from  the  ac  line  is  shaped  into  a  pure  sinusoidal  waveform  with  power  factor 
P.F.=1.0,  problems  such  as  harmonic  currents  and  EMI  will  be  greatly  reduced.  Hence, 
the  power  quality  of  the  utility  source  can  be  improved. 

In  the  whole  analysis  and  design  process,  we  obtained  the  following  new 
contributions.  For  the  active  power  filter  design  in  Chapter  2,  a  new  control  method 
using  an  integrated  and  sampling  technique  is  proposed  to  simplify  the  calculation 
algorithm  for  the  real  fundamental  component  of  load  current.  In  addition,  a  new  simple 
control  scheme  based  on  the  "energy  balance  concept"  is  proposed  to  simplify  the  design 
of  the  conventional  APF-capacitor-voltage-control  circuit.  This  APF  has  many  merits. 
First,  it  simplifies  the  calculation  of  the  current  is(t)  required  by  the  utility  source.  Next, 
a  larger  voltage  ripple  can  be  tolerated  in  the  energy  storage  capacitor;  therefore,  a 
smaller  energy  storage  capacitor  is  needed.  Finally,  only  a  proportional  capacitor- 
voltage-control  method  is  applied  in  the  controller  design,  and  the  transient  response  is 
fast  and  good  (less  than  two  cycles  of  60  Hz). 

In  Chapter  3,  a  new  utility  source  current  memory  circuit  based  on  the  energy- 
balance  concept  is  proposed  to  simplify  the  design  of  the  energy-storage  capacitor  voltage 
control  circuit.  This  design  further  simplifies  the  control  circuit  of  the  "regenerative  load 
system."  For  simplicity,  we  assume  in  Chapter  3  that  the  utility  source  inductance,  Ls, 
equals  zero.  The  effect  of  the  utility  source  inductance  Ls  on  the  system  stability  will  be 
reconsidered  in  Chapter  4  and  is  a  contribution  of  this  research. 


In  the  stability  analysis  of  the  regenerative  load  system,  the  research  in  Chapter  4 
indicated  that  with  the  existence  of  utility  source  impedance  Ls  the  system  will  become 
unstable.  After  the  addition  of  damping  circuit,  the  system  can  be  stabilized;  however, 
the  current  gain  of  the  current  controller  of  the  APF  or  inverter  with  active  power  filtering 
function  must  be  kept  below  one,  i.e.,  k<1.0,  for  the  stability  of  the  system.  If  the  current 
gain,  k,  is  greater  than  one,  then  the  APF  or  inverter  with  active  power  filtering  function 
will  become  a  positive  feedback  system  with  loop-gain  greater  than  one  and  is  an  unstable 
system.  This  phenomenon  has  not  been  noted  in  the  literature. 

Additionally,  in  Chapter  5,  we  derived  the  "input  impedance  specification"  of  the 
DUT  for  the  stability  of  the  regenerative  load  system.  There  are  many  related  research 
works  in  the  literature  [57-68].  However,  we  proposed  an  "equivalent  circuit  concept"  to 
transform  the  system  into  two  subsystems,  an  utility  source  subsystem  and  a  DUT  load 
subsystem.  Then,  the  "two  subsystem  interaction  theory"  by  Schulz  [57]  and  Wildrick  et 
al.  [66]  is  applied  to  analyze  the  stability  of  the  integrated  system.  Based  on  PSpice 
simulation,  we  use  a  simple  current-controller  model  which  is  valid  up  to  1.1  kHz  for  20 
kHz  switching  frequency.  If  the  switching  frequency  is  not  high,  more  detailed  model 
needs  to  be  considered  [70].  This  will  remain  as  an  interesting  future  research  work. 
Finally,  the  input  impedance  specification  of  the  DUT  is  established  for  the  stability  of 
the  regenerative  load  system,  assuming  the  output  impedance  of  the  utility  source 
subsystem  is  known. 

We  also  established  the  "two  damping  circuits  concept."  In  the  analysis  of 
Chapter  5,  we  found  that  if  only  one  damping  circuit  (low-to-medium  frequency)  is  used, 
the  stability  margin  of  the  system  is  not  large  enough.  The  addition  of  a  second  high 
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frequency  damping  circuit  can  further  improve  the  stability  margin  of  the  system.  This 
two  damping  circuit  scheme  can  be  applied  to  other  circuit  designs. 

Lastly,  we  used  the  active  damping  circuit  to  replace  the  passive  damping  circuit 
since  the  passive  damping  circuit  will  dissipate  energy  in  its  resistor  and  decrease  the 
total  efficiency  of  the  system.  Also  these  damping  circuits  can  absorb  the  capacitive- 
reactive  current  and  cause  the  power  factor  of  the  utility  source  to  be  less  than  one,  i.e., 
p.f.  <  1.0.  Therefore,  in  Chapter  6,  the  active  damping  circuit  is  introduced  to  replace  the 
conventional  passive  damping  circuit.  The  energy  consumption  in  the  active  damping 
resistor  is  restored  into  the  energy-storage  capacitor.  Hence  the  efficiency  of  the  system 
can  be  increased.  Also,  in  Chapter  6,  a  "reactive  current  compensation  concept"  is 
developed  to  compensate  the  reactive  current  component  absorbed  by  the  damping 
circuits.  Thus  the  power  factor  of  the  utility  source  current  can  be  increased  to  nearly 
one,  i.e.,  p.f.«  1.0.  In  addition,  it  is  noted  that  due  to  the  bandwidth  limitation  of  the 
current  controller,  only  the  low-to-medium  frequency  damping  circuit  can  be  replaced  by 
the  active  damping  circuit,  and  the  high  frequency  damping  circuit  still  needs  to  be 
passive. 


CHAPTER  2 

ACTIVE  POWER  FILTER 
WITH  REDUCED  ENERGY-STORAGE  CAPACITANCE 

2.1  Introduction 

Power  electronics  are  widely  used  in  industrial  equipment,  such  as  frequency 
changers,  motor-drive  systems,  etc.  Such  equipment  presents  a  nonlinear  impedance 
to  the  utility,  generating  large  harmonic  currents  with  adverse  effects,  such  as  low 
power  factor,  low  efficiency,  and  destruction  of  other  equipment  (e.g.  the  power 
capacitor  can  be  damaged  by  the  resonant  overvoltage,  etc.).  Also,  the  EMI  will 
interfere  with  some  precision  instruments  and  communication  equipment.  Therefore, 
the  active  power  filter  was  used  to  reduce  the  above  pollution  problems  and  to 
compensate  the  reactive  power  of  the  load  simultaneously. 

As  a  result,  utility  power  quality  has  become  an  important  issue.  Many 
research  papers  and  methods  have  been  proposed  to  solve  these  problems. 
Conventionally,  a  passive  LC  filter  was  used  to  suppress  the  harmonics  and  capacitors 
were  used  to  compensate  the  lagging  power  factor.  However,  they  have  many 
disadvantages,  such  as  large  size,  resonance  and  fixed  compensation  characteristics 
[27,28].  The  conventional  passive  power  filter  cannot,  therefore,  provide  a  complete 
solution. 
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Many  specialists  approach  the  solution  from  the  viewpoint  of  preventing  the 
generation  of  harmonics,  such  as  high-power- factor  switching  power  supply,  frequency 
changers  and  uninterruptible  power  supplies  [29-34].  However,  existing  harmonic 
pollution  sources  still  need  to  be  improved.  Some  active  power  filter  (APF)  methods 
have  been  proposed  to  compensate  the  present  harmonic  loads,  to  provide  the  reactive 
power,  and  to  compensate  the  harmonic  current,  in  order  to  achieve  the  goal  of  high 
quality  utility  sources  [25-41].  Since  the  active  power  filter  must  provide  the 
harmonic  and  reactive  power,  it  is  basically  a  bilateral  converter. 

This  chapter  proposes  a  new  APF  technique,  employing  a  reference  voltage  in 
phase  with  the  utility  supply,  to  multiply  with  the  load  current  and  to  calculate  the  real 
fundamental  component  of  the  load  current.  The  difference  between  the  instantaneous 
load  current  and  this  fundamental  component  is  the  command  current  for  the  APF. 
This  command  current  includes  reactive  fundamental  power  and  harmonics.  This 
compensation  current  is  provided  by  a  bilateral  converter  using  the  unipolar  pulse- 
width-modulation  (PWM)  technique  [30,42]. 

As  the  reactive  power  comes  from  the  energy-storage  capacitor  and  this 
reactive  energy  transfers  between  the  load  and  the  energy-storage  capacitor  (charging 
and  discharging  the  capacitor),  the  average  capacitor  voltage  can  be  maintained  at  a 
constant  value.  Owing  to  the  losses  in  the  converter  such  as  switching  loss,  capacitor 
leakage  current,  etc.,  the  utility  must  provide,  not  only  the  real  power  needed  by  the 
load,  but  also  the  additional  power  required  by  the  converter  to  maintain  the  capacitor 
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voltage  at  a  prescribed  value.  In  this  chapter,  we  propose  a  new  dc  bus  voltage 
controller,  using  the  energy  balance  concept. 

As  no  delay  element  (such  as  LPF  or  PI  control)  is  used  in  this  control  circuit, 
the  transient  response  is  fast  and  good  (less  than  two  cycles  of  60  Hz).  The  advantages 
of  this  APF  circuit  are  simplicity  of  control  circuit,  low  cost  (a  smaller  energy-storage 
capacitor),  and  good  transient  response  (less  than  two  cycles  of  60  Hz).  In  theory,  the 
time  delay  for  the  compensation  of  reactive  and  harmonic  currents  is  zero.  The 
feasibility  of  this  theory  is  verified  by  using  a  PSpice  simulation  and  experimental 
results. 

2.2  Principles  of  Operation 

2.2.1.  System  Description 

The  fundamental  building  block  of  the  APF  system  is  shown  in  Fig.  2.1. 
Under  normal  circumstances,  the  utility  can  be  assumed  to  be  a  sinusoidal  voltage 
source. 

vs(t)  =  Vsmsincot  (2.1) 
If  a  nonlinear  load  is  applied,  the  load  current  consists  of  the  fundamental 
component  and  all  the  higher-order  harmonics  and  can  be  represented  as 


iL(t)  =  £lnsin(ncot  +  On) 


(2.2) 


Ufl  i,(t) 


AC  utility 
Vmsin(oi)t) 


mi 


Fig.  2.1  Simplified  schematic  of  the  active-power-line- 
conditioner  or  the  active  power  filter. 
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Therefore  ,  the  load  power  can  be  expressed  as 


PL(t)  =  vs(t)iL(t) 

=  IiVsmsin  cotcosG] 
+I|Vwlsiil  cotcosG! 

oo 

+2  Vsm  sincotIn  sin(ncot  +  9n) 

n=2 

=  Ps(t)  +  Pc(t)  (2.3) 
In  Eq.  (2.3),  the  first  term  is  the  real  power  supplied  by  the  utility  source,  the  second 
term  is  the  reactive  power,  and  the  third  term  is  the  harmonic  power,  i.e., 


Ps(t)  =  IiVsm  sin2  cotcos0, 


pc(t)  =  1 1  Vsm  sincot  cos  cot  sin  0! 


oo 


+  ]TVsmsin(cot)Insin(ncot  +  en)  (2.4) 

n=2 

As  the  APF  provides  the  reactive  and  harmonic  power  pc  (t) ,  the  current  supplied  by 
utility  will  be 

h(t)  =        =  I,  cosG,  sincot 
vs(t) 


=  IS  sincot 


where 


I,  =  I,  cosB,  (2.5) 


The  current  is  in  phase  with  the  utility  voltage  and  is  pure  sinusoidal.  At  this  time,  the 
APF  must  provide  the  following  compensation  current 


ic(t)  =  iL(t)-is(t)  (2.6) 

Thus  the  APF  needs  to  calculate  the  current  is(t)  accurately  and  instantaneously  as 
described  below. 

2.2.2  Calculation  of  the  Fundamental  Component  u(t) 

Now,  if  a  sinusoidal  waveform  with  amplitude  2  and  in  phase  with  the  voltage 
source  vs(t)  is  multiplied  by  iL(t), 

00 

p'L(t)  =  2  sin  cot  J]  I  n  sin(ncot  +  9n)  ....  (2.7) 

n=l 

Integrating  Eq.  (2.7)  over  a  period  T,  the  averaging  power  can  be  obtained  as 


Pav=:j:(pL(t)dt 


(21 1  sin2 cotcos9,)dcot 


CO  J 

+^ —  |(2sincotIn  sin(ncot  +  Gn)dcot 

n=2  frt 


~coT 

1  |<oT 

+  — I    (21,  sincotcoscotsinSi  )dcot 
coT  * 

a>T 

IiCosGj  (2.8) 
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After  integration  over  a  whole  period,  the  second  and  third  terms  on  the  right-hand 
side  of  Eq.  (2.8)  will  be  equal  to  zero.  Therefore,  Pav  =  I,  cosGj  is  the  peak  value  of 
the  current  supplied  by  the  voltage  source  ,  as  seen  in  Eq.  (2.5). 

The  block  diagram  of  the  fast  utility  source  current  calculation  circuit  is  shown 
in  Fig.  2.2.  The  voltage  source  is  buffered  down  and  fed  through  a  zero-cross  detector 
(ZCD)  to  provide  a  synchronous  square  wave  as  the  input  of  the  phase-locked-loop 
(PLL)  circuit.  The  PLL  circuit  will  produce  a  reference  signal  2sincot  in  phase  with  the 
utility  voltage  source.  This  signal  is  multiplied  by  the  load  current  i  L  (t) ,  integrated 
over  a  period  T,  and  fed  through  a  sample-and-hold  (S/H)  circuit  to  obtain  a  quantity 
proportional  to  Ij  cos^ .  The  control  logic  resets  the  integrator  at  the  beginning  of 

every  period  T  and  samples  the  integration  results  at  the  end  of  the  integration.  In  this 
way,  the  real  power  can  be  calculated.  By  appropriate  adjustment  of  the  gain  K,  the 
peak  value  of  the  utility  source  current  can  be  obtained. 

2.2.3.  Control  of  the  dc-bus  Capacitor  Voltage 

The  compensation  current  ic(t)  is  provided  by  a  bilateral  converter  with 
unipolar  PWM  modulation  as  shown  in  Fig.  2.3.  The  dc  source  comes  from  the  energy 
storage  capacitor,  as  it  provides  only  the  reactive  power  to  the  load.  However, 
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Vs(0 


Differential 
amplifier 


Zero  cross 
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phase-locked 
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Control  logic 
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Integrator 


sampling 
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Voc(t) 
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Fig.  2.2  Block  diagram  of  the  fast  calculation  circuit  for  the  utility  source  current. 
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Fig.  2.3  The  bilateral  converter  using  the  unipolar  PWM  modulation  to  provide  the 
compensation  current  ic(t). 
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owing  to  the  switching  loss  of  the  converter,  the  utility  must  supply,  not  only  the  real 
power  of  the  load  current,  but  also  a  small  overhead  for  the  capacitor  leakage  and 
converter  switching  losses.  In  this  way,  the  dc-bus  voltage  can  be  maintained  on  a 
fixed  level.  The  control  strategy,  using  the  energy-balance  principle,  is  described  in 
the  following. 

If  the  reference  voltage  across  the  capacitor  is  Vcr ,  the  reference  energy  storage 

will  be 

Ecr  =  ^CcVc2r  (2.9) 
while  the  instantaneous  energy  in  the  capacitor  is 

Ec(t)  =  ^Ccvc2(t)  (2.10) 
Therefore,  the  energy  loss  of  the  capacitor  in  one  period  will  be 

AEc(t)  =  Ecr-Ec(t) 

=  ^{Vc2r-vc2(t)} 

=  ^{Vcr+vc(t)}{Vcr-vv(t)}  (2.11) 

Assuming  that,  within  one  period,  the  variation  in  capacitor  voltage  is  moderate,  the 
term  {Vcr  +  vc(t)}  can  be  approximated  as 
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Vcr+vc(t)*2Vcr  (2.12) 

AEc(t)*CcVcr{Vcr-vc(t)} 

=  K,{Vcr-vc(t)}  (2.13) 
Since  this  energy  loss  must  be  supplied  by  the  utility  voltage  source,  the  peak  value  of 
the  charging  current  Isc  can  be  estimated  as 


f  Vc 


sm 


Therefore 


ISC=-|-AEC=K2AEC 

1  V  c 


sm 


where 

K2  = 


TV 

1  vsm 


Hence,  the  peak  value  of  the  charging  current  Isc  can  be  obtained  as 


Isc=K2AEc  =  K2K,{Vcr-vc(t)} 
=  K!2{Vcr-vc(t)} 


where 


(2.14) 


Ki2=K,K2=^S-  (2.15) 
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As  a  result,  the  total  peak  current  provided  by  the  voltage  source  is 


Ism=Is+I 


sc 


is(t)  =  Ism  sincot 


(2.16) 


C(t)  =  iL(t)-is(t) 


(2.17) 


By  multiplying  Ism  with  sincot,  we  can  obtain  the  utility  source  current  is(t).  The 
difference  between  the  load  current  iL(t)  and  is(t)  will  be  the  command  current 


Either  bipolar  or  unipolar  voltage  switching  of  the  converter  can  be  used.  The 
latter  is  selected  to  reduce  converter  loss.  Thus,  the  current  tracking  circuit  is  a 
bilateral  converter  with  unipolar  PWM  technique  [30,42]  as  shown  in  Fig.  2.3.  The 
switching  frequency  is  25kHz. 


Figure  2.4  shows  the  simulation  results  using  PSpice  software,  where 
Li=2.5mH,  Ci=10|iF,  Cc=660u.F,  and  f s  =  25kHz.  When  the  voltage  source  is 
applied,  the  estimation  integration  circuit  for  the  current  is(t)  has  not  yet  provided  the 
calculated  value.  Hence,  the  output  is  zero,  and  all  the  load  current  is  supplied  by  the 
converter.  The  capacitor  voltage  will  drop  owing  to  the  release  of  stored  energy.  In 
the  second  cycle,  however,  the  utility  will  supplement  the  energy  in  the  capacitor  to 


i*  (t)  for  the  APF  converter. 


2.3  Simulation  Results 
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Fig.  2.4  Simulated  current  and  voltage  waveforms  for  steady-state  and  transient  response. 

(a)  Utility  source  voltage 

(b)  Load  current 

(c)  Utility  source  current 

(d)  Compensation  current 

(e)  Capacitor  (Cdc)  voltage 

Load  current  changes  from  a  peak  value  of  10A  to  20A  after  the  fifth  cycle. 
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raise  its  voltage  to  the  prescribed  value.  Even  so,  the  utility  current  still  maintains  a 
pure  sinusoidal  waveform.  Except  for  the  one-period  delay  by  the  integrator  (this  is 
the  minimum  time  required  to  calculate  the  real  power),  the  whole  control  circuit  has 
no  other  time  delay.  Therefore,  the  transient  response  will  be  fast  and  good.  In  the 
simulation  circuit,  there  is  a  step  change  of  load  current,  from  a  peak  value  of  10A  to 
20A,  after  the  fifth  cycle.  However,  the  utility  source  current  responds  immediately  at 
the  next  cycle  and  stabilizes  at  the  new  value.  This  behavior  demonstrates  the  fast  and 
good  transient  response. 

Figure  2.5  shows  the  frequency  components  of  load  current  iL  utility  source 
current  is,  and  converter  compensation  current  ic  obtained  by  PSpice  simulation.  In 
Fig.  2.6,  the  transient  response  of  a  step  change  of  load  current  from  peak  value  of 
20 A  to  10A  at  the  fifth  cycle  is  presented.  From  this  simulation,  a  capacitor- voltage 
rise  during  the  transition  cycle  can  be  observed. 

If  limitation  of  the  capacitor  voltage  is  required,  the  utility  current  must  change 
just  after  the  voltage  boundary  is  encountered.  To  achieve  this  an  additional  circuit 
can  be  added  to  detect  the  crossing  of  the  voltage  limit  and  to  estimate  the  extra 
fundamental  utility  current  needed  to  return  the  capacitor  voltage  to  the  preset  value 
before  the  beginning  of  the  next  cycle. 

Figure  2.7  shows  the  results  of  PSpice  simulation,  where  the  voltage  limit  is  set 
as  vmax  =  235V  and  Vmin  =  175V .  In  addition,  the  load  current  changes  from  peak 
value  of  10A  to  20A  at  the  fifth  cycle  and  from  peak  value  of  20A  to  10A  at  the  tenth 
cycle. 
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Fig.  2.5  Frequency  components  of  load  current  iL ,  utility  source  current  is  and 
converter  compensation  current  ic 

(a)  Load  current  i  L  =  i(vm) 

(b)  Utility  source  current  is=i(vms) 

(c)  converter  compensation  current  ic=i(vmc) 
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Fig.  2.6  Simulated  current  and  voltage  waveforms  for  steady-state  and  transient  response. 

(a)  Utility  source  current 

(b)  Load  current 

(c)  Utility  source  current 

(d)  Conpensation  current 

(e)  Capacitor  (Cdc)  voltage 

Load  current  changes  from  a  peak  value  of  20A  to  10A  after  the  fifth  cycle. 
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Fig.  2.7  Simulated  current  and  voltage  waveforms  for  steady-state  and  transient  response. 

(a)  Utility  source  voltage 

(b)  Load  current 

(c)  Utility  source  current 

(d)  Compensation  current 

(e)  Capacitor  (Cdc)  voltage 

Maximum  and  minimum  voltages  across  the  capacitor  are  preset  at  Vmax  =  235V  and 
Vmin=  175  V 
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2.4  Experimental  Results 

Figure  2.8(a)  shows  the  experiment  results  where  the  RMS  value  of  voltage 
source  is  28v  (voltage  waveform  pre-scaled  down  by  1/3),  the  load  on  the  dc  side  is  a 
4700uF  capacitor  in  parallel  with  36Q  resistor,  and  idc  =  0.85A ,  whereas  Fig.  2.8(b)  is 
the  result  when  a  40uF  capacitor  is  placed  in  front  and  parallel  with  the  rectifier.  The 
dc-bus  voltage  of  the  converter  is  50v  and  Cc  is  330uF. 

Figure  2.9(a)  and  (b)  are  the  experiment  results  for  transient  response.  The 
load  current  began  with  a  value  of  1  p.u.  RMS,  rose  up  to  2  p.u.  and  fell  back  again  to 
1  p.u.  In  Fig.  2.9(a),  the  energy-storage  capacitor  is  assumed  as  1.42  p.u.,  whereas  in 
Fig.  2.9(b),  the  capacitor  value  is  0.46  p.u.,  and  the  voltage  across  capacitor  "Cc"  is  set 
at  1.76  p.u.  From  Fig.  2.9(b),  it  can  be  seen  that  even  with  a  large  ripple  voltage  on 
the  energy-storage  capacitor,  the  compensation  results  are  still  good. 

2.5  Discussion 

The  determination  of  the  value  of  the  energy-storage  capacitor  is  based  on  the 
following  three  conditions: 

(1)  Step  increase  of  the  real  fundamental  component  of  the  load  current 

When  the  load  current  has  a  step  increase,  the  energy  stored  in  the  capacitor 
must  be  released  immediately  to  support  the  step  increase  in  the  power  consumed  by 
the  load,  i.e.,  by  using  energy-balance  concept 
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Fig.  2.8  Experimental  results  for  steady-state  where  Cc=  330uF:  (a)  Load  on  dc  side 

is  a  4700uF  capacitor  in  parallel  with  a  36  Q  resistor;  (b)  40  uP  capacitor  in  front  and 
parallel  with  the  rectifier 

(1)  Utility  voltage  (voltage  waveform  pre-scaled  down  by  1/3) 

(2)  Load  current 

(3)  Utility  source  current 

(4)  Compensation  current 
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(b) 


Fig.  2.9  Experimental  transient-response  results:  (a)  Cdc=1.42  p.u.;  (b)  Cdc=0.46  p.u. 

Ccbasc  =1P-u-=(IbaseT)/Vbasc 

Vbase  =  Vs.RMS 
'base  =  h.RMS 

T=60Hz 

(1)  Utility  voltage  (voltage  waveform  pre-scaled  down  by  1/3) 

(2)  Load  current 

(3)  Utility  source  current 

(4)  Voltage  across  energy-storage  capacitor  Cc 
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^Cc{Vc2r-Vc2min}=^VsmAlL1T  '  (2.18) 


where  Vc  min  =  preset  lower  limit  of  the  voltage  of  the  energy-storage  capacitor,  AIL, 

=  step  increase  of  the  peak  value  of  the  real  fundamental  component  of  the  load 
current  and  T  =  period  of  the  utility  voltage  source.  Therefore,  the  value  of  the 
capacitor  is 


^ai^j_  (219) 

v  -  V 

Tcr  Tc,min 

(2)  Step  reduction  of  the  real  fundamental  component  of  the  load  current 

While  the  load  current  is  reduced,  the  utility  source  current  does  not  change 
until  the  next  cycle.  Hence,  this  extra  utility  source  current  AIL2  will  charge  the 
energy-storage  capacitor,  i.e.,  by  using  the  energy-balance  concept 


^Cc{Vc2max-Vc2r}  =  ivsmAIL2T  (2.20) 

where  Vc  max  =  preset  upper  limit  of  the  voltage  of  the  energy  storage  capacitor  and 

AI L2  =  step  reduction  of  the  peak  value  of  the  real  fundamental  component  of  load 
current.  Therefore,  the  value  of  the  capacitor  required  is 
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(3)  Reactive  and  harmonic  components  of  the  load  current 

During  the  steady  state,  the  reactive  and  harmonic  components  of  the  load 
current  will  charge  and  discharge  the  energy-storage  capacitor  during  the  period. 
Therefore,  using  the  energy  balance  concept 


2  c 


(Vc2A-Vc2r)=^VsmAIL3^-  (2.22) 


1-.  T 


where  AI L3  =  the  peak  value  of  the  RMS  of  the  reactive  and  harmonic  component  of 
load  current  and  VcA  =  maximum  or  minimum  voltage  of  the  energy-storage  capacitor 
during  one  period.  Therefore,  using  the  energy  balance  concept 


VsmAIL3T 


Cc3  = 


(vc2A-vc2r) 


(2.23) 


where  the  half  period  T/2  is  selected  as  the  maximum  charging  or  discharging  times  of 
the  reactive  and  harmonic  components  will  be  less  than  T/2. 

If  the  sinusoidal  current  waveform  during  the  transient  is  to  be  preserved,  the 
largest  of  Ccl  ,  Cc2  ,  and  Cc3  is  selected  as  the  value  of  the  energy-storage  capacitor. 
If  the  sinusoidal  current  waveform  is  not  needed  during  the  transient,  the  value  of  Cc3 
is  chosen  as  the  value  of  the  energy-storage  capacitor.  In  Fig.  2.7,  the  latter  case  is 
applied  in  order  to  reduce  the  value  of  the  maximum  voltage  across  the  energy-storage 
capacitor. 
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2.6  Summary 

Owing  to  the  nonlinear  load  of  much  electronic  equipment,  the  utility  power 
system  is  polluted  by  harmonics  and  the  power  factor  is  reduced.  In  this  chapter,  a 
new  APF  technique  is  proposed,  using  integration  of  the  power  and  a  sampling 
technique  to  simplify  the  calculation  of  the  real  fundamental  component  of  load 
current.  In  addition,  the  energy-balance  concept  in  the  energy-storage  capacitor  is  used 
to  simplify  the  design  of  the  conventional  APF-capacitor-voltage-control  circuit.  This 
APF  has  three  merits.  To  begin  with,  the  calculation  circuit  for  the  utility  source 
current  is(t)  is  simplified.  Second,  a  larger  ripple  voltage  can  be  tolerated  in  the 
energy  storage  capacitor;  therefore,  a  smaller  energy  storage  capacitor  is  needed. 
Third,  only  a  proportional  capacitor-voltage-control  method  is  applied  in  the  controller 
design.  Hence,  the  transient  response  is  fast  and  good  (less  than  two  cycles).  The 
feasibility  of  the  above  scheme  is  verified  by  PSpice  simulation  and  experimental 
results.  The  results  demonstrate  that  harmonics  are  suppressed  and  a  nearly  unity 
power  factor  is  obtained. 


CHAPTER  3 
REGENERATIVE  LOAD  SYSTEM 

3.1  Introduction 

For  testing  the  power  device,  if  a  conventional  load  is  used,  a  large  amount  of 
energy  will  be  dissipated  and  Joule  heat  will  be  produced.  Hence,  the  focus  of  this 
chapter  is  trying  to  develop  a  regenerative  load  for  long  term  testing  of  power 
equipment  and  motor-generator  set  to  achieve  energy  conservation  and  improve  the 
power  quality.  The  proposed  system  is  shown  in  Fig.  3.1.  Basically  this  system  will 
consist  of  two  parts:  an  impedance  emulator  and  an  inverter  with  active  power  filtering 
function.  The  internal  structural  block  diagram  is  shown  in  Fig.  3.2. 

In  the  first  place,  the  impedance  emulator  draws  current  from  the  DUT, 
according  to  the  equation  Il=Vl/Il.  This  current  related  power  can  be  used  to  charge 
the  energy  storage  capacitor;  therefore,  the  input  of  this  converter  can  be  viewed  as  an 
equivalent  load  impedance  of  value  Zl 

The  energy  stored  in  the  capacitor  is  inverted  into  ac  current  and  fed  back  into 
the  input  of  the  DUT  by  the  second  part  of  the  system,  i.e.,  the  inverter  with  active 
power  filtering  function.  In  addition,  this  inverter  can  also  provide  the  compensation 
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Fig.  3.1  The  regenerative  load  system. 
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Fig.  3.2  The  internal  structural  block  diagram  of  the  regenerative  load 
system 
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current  to  shape  the  ac  line  current  is  into  pure  sinusoidal  waveform  with  power  factor 
P.F.=1.0. 

Hence,  the  proposed  system  has  three  merits.  Because  the  system  recycles  the 
electrical  energy,  the  energy  consumption  will  be  less  than  20%  of  the  conventional 
load.  If  the  soft-switching  scheme  is  employed  for  the  converters,  the  dissipation  can 
be  even  lower  than  10%.  Additionally,  the  system  is  safe.  Due  to  the  low  energy 
dissipation,  no  appreciable  Joule  heating  is  observed.  Lastly,  since  the  current  from 
the  ac  line  is  shaped  into  sinusoidal  waveform  with  power  factor  P.F.=1.0,  problems 
such  as  harmonic  currents  and  EMI  will  be  greatly  reduced.  Thus,  the  power  quality 
of  the  utility  source  can  be  improved. 

In  this  chapter,  a  new  dc  capacitor  voltage  controller  using  a  new  utility  source 
current  memory  circuit  and  the  energy-balance  concept  is  proposed.  In  order  to 
concentrate  on  the  description  of  the  energy-balance  concept  and  the  utility  source 
current  memory  circuit,  we  assume,  temporarily,  the  utility  source  is  an  ideal  voltage 
source,  i.e.,  the  source  inductance,  Ls,  is  zero.  The  effect  of  utility  source  inductance, 
Ls,  on  the  system  stability  will  be  discussed  in  Chapter  4. 

3.2  Principles  of  Operation 
The  example  schematic  of  the  regenerative  load  system  is  shown  in  Fig.  3.3, 
and  consists  of  the  impedance  emulator  and  the  inverter  with  active  power  filtering 
function. 


34 


Fig.  3.3  The  schematic  circuit  diagram  of  the  regenerative  load  system. 


3.2.1  Impedance  Emulator 

For  the  impedance  emulator,  a  simple  dc-to-dc  boost  converter  is  employed  as 
a  dc  current  controller  with  the  reference  current  such  that 

I*ef  =  V  o  /  Rl  (3.1) 

where  Vo  is  the  output  voltage  of  the  device  under  test  and  Rl  is  the  emulated  value  of 
the  load.  Therefore,  the  output  current  of  the  device  under  test,  Io,  will  track  the 
reference  current,  I  ref .  The  energy  into  the  impedance  emulator  is  used  to  charge  the 

capacitor  Cdc.  The  control  circuit  of  the  impedance  emulator  is  shown  in  Fig.  3.4. 
Switches  Seland  Se2are  used  to  select  a  reference  current  I  ref .  Once  the  reference 
current  Iref  is  selected,  the  actual  current  into  the  boost  converter  will  track  the 
reference  value  I  ref .  The  current  tracking  scheme  used  is  conventional  bang-bang 
hysteresis  controller.  The  output  signal  is  used  to  control  the  power  stage  switch  Sie 
of  the  impedance  emulator  shown  in  Fig.  3.3. 

3.2.2  Invertor  with  Active  Power  Filtering  Function 

Calculation  of  the  command  current  ig*(t)  for  the  inverter 

As  for  the  second  part,  the  inverter  with  active  power  filtering  function  consists 
of  a  single-phase  bridge  converter  acting  as  a  single-phase  current  controller.  The 
voltage  across  the  capacitor  Cdc  acts  as  the  voltage  source  of  this  inverter. 


Fig.  3.4  The  control  circuit  of  the  impedance  emulator. 


37 

The  command  current  ic*(t)  for  this  inverter  is  calculated  as  shown  in  Fig.  3.5. 
The  control  strategy,  using  the  energy-balance  principle,  is  described  below. 

If  the  reference  voltage  across  the  capacitor  is  Vcr,  the  reference  energy  storage 

will  be 


Ecr  =  ICcVc2r  (3.2) 


while  the  instantaneous  energy  in  the  capacitor  is 


1 


Ec(t)  =  ^Ccvc2(t)  (3.3) 


2 

Therefore,  the  energy  loss  of  the  capacitor  in  one  period  will  be 


AEc(t)  =  Ecr  -  Ec(t) 
Cc 


_c  r-w2  ..2, 

=  y{Vcr  +  Vc(t)}{Vcr-Vc(t)}  (3.4) 

Assuming  that,  within  one  period,  the  variation  in  capacitor  voltage  is  moderate,  the 
term  {Vcr  +  ve(t)}  can  be  approximated  as 


Vcr  +  Vc(t)  «  2Vcr  (3.5) 

AEc(t)  *  CcVcr{Vcr  -  Vc(t)} 

=  Kl{Vcr- Vc(t)}  (3.6) 
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Fig.  3.5  The  block  diagram  of  the  calculating  circuit  for  the  utility  source  current. 


Since  this  energy  loss  must  be  supplied  by  the  utility  voltage  source,  the  change  of  the 
peak  value  of  the  charging  current,  A  Isc,  can  be  estimated  as 


f 


Vsm  sincotAIsc  sincotdt  =  AEc 


(3.7) 


Therefore 


Alsc  =  AEc  =  K.2  AEc 

TVsm 


(3.8) 


where 


K2  = 


TVs 


Hence,  the  variation  of  the  peak  value  of  the  charging  current,  Alsc,  can  be  obtained  as 
follows 


Alsc  =  K2AEc  =  K2Kl{Vcr  -  vc(t)} 


=  Kl2{Vcr-Vc(t)} 


(3.9) 


where 


Ki2=  K2Ki  = 


2CcVcr 


TV 


sm 


In  the  next  period,  if  this  amount  of  peak  value  A  Isc  =  A  Ii  is  added  to  the  steady-state 
utility  current  Isci,  the  capacitor  voltage  will  be  restored  by  an  amount  of  A  V|  as 
shown  in  Fig.  3.6.  However,  the  steady-state  utility  source  current  must  also  be 
updated  by  an  amount  AI2  =  Alsc  in  order  to  prevent  capacitor  voltage  from  continuing 
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(a)  Voltage  change  of  the  energy  storage  capacitor  Cdc. 


*sc2 


|ai, 

\**2 


lsc3 


-*■  t 


\  / 

sampling  at  zero 
crossing  of  60Hz 


(b)  Peak  value  of  the  charging  current  I 


Fig.  3.6  The  concept  for  the  memory  of  the  peak  value  of  steady-state  utility  source 
current,  ig. 

1 .  AIj  is  used  to  adjust  the  steady-state  charging  current  Isc  in  order  to  prevent 
vc(t)  from  continuing  change  (e.g.,  voltage  drop  AV2). 

2.  AI,  is  used  to  raise  the  vc(t)  for  an  amount  AV,  =  AV. 

3.  AV,  =  AV2  =  AV  and  AI,  =  AIj  =  AIK. 

4.  1^,  is  the  previous  memory  of  the  peak  value  of  steady-state  utility  current. 

5.  1^  is  the  updated  memory  of  the  peak  value  of  steady-state  utility  current. 

6.  is  the  actual  peak  value  of  1^  during  t,  ~  tj. 
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to  change  (e.g.  another  voltage  drop  AV2 ).  Therefore,  the  memory  of  the  peak  value 
of  the  updated  steady-state  utility  current  needs  to  be 

ISC3  =Iscl  +AI2  (3.10) 

where 

Al2  =  AIscSKl2{Vcr-Vc(t)} 

As  a  result,  during  the  next  period,  the  total  peak  value  of  source  current  will  be 


Isc2=AI,+Isc3 


=  AI,  +AI2  +ISC| 

=  2AI  +  Iscl  (3.11) 


where 


AI  =  All  =  AI2  =  Kl2{Vcr  -  Vc(t)} 

By  multiplying  Isc2  with  sincot,  it  is  possible  to  obtain  utility  source  current  is(t).  The 
difference  between  the  load  current  it(t)  and  is(t)  will  be  the  command  current  ic('t)  for 
the  inverter  with  active  power  filtering  function. 


is(t)  =  Isc2  sincot  (3.12) 


ic(t)  =  iL(t)-is(t)  (3.13) 

During  steady-state,  the  capacitor  voltage  vc(t)  at  the  sampling  instant  will  equal  to 
reference  voltage  Vcr. 
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Vc  (0  sampling  =Vcr  (3.14) 

Hence,  the  variation  of  the  charging  current  AI  will  be 

AI  =  AIsc  =  Kl2{Vcr-Vc(t)}  =0  (3.15) 
Consequently,  in  the  steady-state,  the  utility  source  current  will  become 

Isc2  =  2AI  +  Iscl  =  Iscl 

is(t)  =  Isc2  sincot  =  Isci  sincot 

=  Iss  sincot  (3.16) 

where 

Iss  =  steady-state  peak  value  of  utility  source  current 

In  addition,  the  updated  memory  of  ISC3,  will  be  re-sampled  as  the  previous  memory 
for  the  following  cycle  as  shown  in  Fig.  3.5. 

The  detailed  schematic  diagram  of  the  calculation  circuit  of  the  command 
current  ic*(t)  for  the  inverter  with  filtering  function  is  shown  in  Fig.  3.7. 
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Current-mode  control  scheme  of  the  inverter 

Either  bipolar  or  unipolar  voltage  switching  of  the  inverter  can  be  used.  The 
latter  is  selected  to  reduce  the  inverter  loss.  Hence,  the  current  tracking  circuit  is  a 
bipolar  converter  with  unipolar  PWM  technique  [30,  42],  as  shown  in  Fig.  3.8.  The 
detailed  control  circuit  is  shown  in  Fig.  3.9,  and  the  switching  frequency  is  25  kHz. 

3.3  Simulation  Results 
Figure  3.10  shows  the  simulation  results,  using  PSpice  software,  where  Li  = 
2.5mH,  Ci  =  lOuF,  Cdc  =  660uF  (or  470uF)  and  fs  =25  kHz.  When  the  utility  source 
current  is  applied,  the  impedance  emulator  has  not  been  turned  on.  Hence,  all  the 
currents  are  nearly  zero  and  the  capacitor  voltage  is  kept  at  the  reference  voltage  Vref 
=  285v.  After  the  fifth  cycle,  the  impedance  emulator  is  turned  on,  and  there  is  a  step 
change  of  the  output  current  of  the  DUT  from  i0  =  OA  to  i0  =  10A.  The  input  current  to 
the  DUT  is  stabilized  at  the  seventh  cycle;  however,  the  capacitor  voltage  vdc  is  lower 
than  Vref .  Hence,  the  capacitor  voltage  is  charged  to  Vref  in  the  eighth  cycle.  The 
whole  system  is  stabilized  at  the  beginning  of  the  ninth  cycle.  Therefore,  the  transient 
response  is  two  cycles  longer  than  that  of  DUT.  When  the  DUT  is  stabilized,  the 
capacitor  voltage  may  not  be  at  the  Vref  value.  As  a  result,  there  is  one  more  cycle 
needed  for  the  capacitor  voltage  to  be  charged  to  Vref . 
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Fig.  3  .8  The  bilateral  current  controller  using  the  unipolar  PWM  modulation. 


Fig.  3  .9  The  detailed  control  circuit  of  the  bilateral  current  controller  with 
unipolar  PWM  modulation. 
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Fig.  3.10  Simulated  current  and  voltage  waveforms  for  the  regenerative  load  system. 

(a)  Utility  source  voltage  vs  =  v(7,3) ;  utility  source  current  is  =  i(vms) 

(b)  DUT  input  current  iL  =  i(vm) 

(c)  Current  controller  current  icc  =  i(vmc) 

(d)  DUT  output  voltage  v0  =  v(l) ;  DUT  output  current  i0  =i(vmle) 

(e)  Capacitor  (Cdc)  voltage  Vdc  =  v(71,FGnd) 

DUT  output  current  changes  from  a  peak  value  of  OA  to  10A  at  the  fifth  cycle. 


In  Fig.  3.1 1,  the  transient  response  of  a  step  change  of  DUT  output  current 
from  i0(t)  OA  to  i0(t)  -  OA  is  presented.  From  this  simulation,  a  two  cycle  longer 
transient  response  than  that  of  DUT  can  be  observed.  Figure  3.12  shows  the  results  of 
PSpice  simulation,  where  the  output  current  of  DUT  changes  from  i0(t)  =  OA  to  10A  at 
the  fifth  cycle  and  from  10A  to  OA  at  the  tenth  cycle. 


The  determination  of  the  value  of  energy-storage  capacitor  is  based  on  the 
following  three  conditions: 

(1)  Step  increase  of  the  real  fundamental  component  of  the  DUT  input  current 

When  the  DUT  input  current  has  a  step  increase,  the  energy  stored  in  the 
capacitor  must  be  released  immediately  to  support  the  step  increase  in  the  power 
consumed  by  the  DUT,  i.e.,  by  using  the  energy-balance  concept 


where  Vc,min  =  preset  lower  limit  of  the  voltage  of  energy-storage  capacitor,  AIli  a 
step  increase  of  the  peak  value  of  the  real  fundamental  component  of  the  DUT  input 
current  and  T  =  period  of  the  utility  voltage  source.  Therefore,  the  value  of  the 
capacitor  is 


3.4  Discussion 


5<*0£  - 


v.2min}  =  -vsmAlL1T 
c,minJ  2 


(3.17) 
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Fig.  3.11  Simulated  current  and  voltage  waveforms  for  the  regenerative  load  system. 

(a)  Utility  source  voltage  vs  =  v(7,3) ;  utility  source  current  is  =  i(vms) 

(b)  DUT  input  current  iL  =  i(vm) 

(c)  Current  controller  current  icc  =  i(vmc) 

(d)  DUT  output  voltage  v0  =  v(l) ;  DUT  output  current  i0  =  i(vmle) 

(e)  Capacitor  (Cdc)  voltage  Vdc  =  v(71  ,FGnd) 

DUT  output  current  changes  from  a  peak  value  of  1  OA  to  OA  at  the  fifth  cycle. 
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Fig.  3.12  Simulated  current  and  voltage  waveforms  for  the  regenerative  load  system. 

(a)  Utility  source  voltage  vs  =  v(7,3) ;  utility  source  current  is  =  i(vms) 

(b)  DUT  input  current  iL  =  i(vm) 

(c)  Current  controller  current  icc  =i(vmc) 

(d)  DUT  output  voltage  v0  =  v(l) ;  DUT  output  current  i0  =i(vmle) 

(e)  Capacitor  (Cdc)  voltage  Vdc  =  v(71,FGnd) 

DUT  output  current  changes  from  a  peak  value  of  OA  to  10A  at  the  fifth  cycle  and  from 
10A  to  OA  at  the  tenth  cycle. 
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(2)  Step  increase  of  the  real  fundamental  component  of  the  DUT  output  current 

While  the  DUT  output  current  is  increased,  its  input  current  may  not  change  as 
quickly  and  as  large  until  the  next  cycle.  Hence,  this  extra  DUT  output  current  AIL2 
will  charge  the  energy  storage  capacitor,  i.e.,  by  using  the  energy  balance  concept 


^Cc{V2max  -  V2}  =  lvsmAIL2T  (3.19) 


where  Vc,max  =  preset  upper  limit  of  the  energy  storage  capacitor  and  AIL2  =  step 
increase  of  the  peak  value  of  the  real  fundamental  component  of  DUT  output  current. 
Therefore,  the  value  of  the  capacitor  required  is 

Cc2=  VsmAlL2T  (32Q) 

V2  -  V2 
Y  cmax     T  cr 

(3)  Reactive  and  harmonic  components  of  the  DUT  input  current 

During  the  steady-state,  the  reactive  and  harmonic  components  of  the  DUT 
input  current  will  charge  and  discharge  the  energy-storage  capacitor  during  the  period. 
Therefore,  using  the  energy  balance  concept 


ICC 
2 


(vc2A-vc2r) 


=  lVsmAIL3|  (3.21) 
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where  AIL3  =  the  peak  value  of  the  RMS  of  the  reactive  and  harmonic  component  of 
DUT  input  current  and  Vca  =  maximum  or  minimum  voltage  of  the  energy-storage 
capacitor  during  one  period.  Therefore,  using  the  energy  balance  concept 


T 

VsmAIL3  — 

Cc3  =  -—  2-  (3.22) 

(VC2A-Vc2r) 

where  the  half  period  172  is  selected,  so  that  the  maximum  charging  or  discharging 
times  of  the  reactive  and  harmonic  components  will  be  less  than  T/2. 

If  the  transient  response  waveform  of  the  system  is  to  be  preserved,  the  largest 
of  Cci,  Cc2  and  Cc3  is  selected  as  the  value  of  the  energy-storage  capacitor.  If  the 
transient  response  waveform  of  the  DUT  is  not  needed,  the  value  of  Cc3  is  chosen  as 
the  value  of  the  energy-storage  capacitor. 


3.5  Summary 

For  the  long  term  testing  of  power  equipment  and  motor-generator  set,  a 
regenerative  load  system  is  proposed,  based  on  the  energy  balance  concept.  In 
addition,  a  memory  circuit  is  added  to  the  dc  capacitor  voltage  controller  to  store  the 
value  of  the  utility  source  current  needed  to  maintain  the  capacitor  voltage  during 
steady-state.  This  memory  circuit  concept  can  simplify  the  design  of  the  capacitor 
voltage  controller.  This  regenerative  load  system  has  three  advantages:  (1)  economic, 
the  energy  consumption  will  be  less  than  20%  of  the  conventional  load;  (2)  safety,  due 


to  the  low  energy  dissipation,  no  appreciable  Joule  heating  is  observed;  and  (3)  power 
quality,  since  the  ac  line  current  is  shaped  into  sinusoidal  waveform  with  P.F.  =  1.0, 
problems  such  as  harmonic  currents  and  EMI  are  greatly  reduced  and  the  power 
quality  of  the  utility  source  is  also  improved.  In  addition,  only  a  proportional 
capacitor-voltage-control  method  is  applied  in  the  controller  design.  Thus,  the 
transient  response  is  fast  and  good.  The  feasibility  of  this  regenerative  load  system  is 
verified  by  PSpice  simulation. 


CHAPTER  4 

STABILITY  ANALYSIS  OF 
REGENERATIVE  LOAD  SYSTEM 
WITH  UTILITY  SOURCE  INDUCTANCE 

4.1  Introduction 

For  the  regenerative  load  system  in  Chapter  3,  we  assumed  that  the  utility 
source  is  an  ideal  voltage  source  and  has  no  series  inductance  with  itself.  However, 
actually,  there  will  be  some  leakage  inductance  of  the  power  transformer,  Ls,  in  series 
with  the  utility  source.  When  this  series  inductance  is  taken  into  consideration,  the 
regenerative  load  system  will  become  marginally  stable  or  unstable. 

This  chapter  will  discuss  the  stability  problem  of  the  regenerative  load  system 
with  the  source  inductance  taken  into  consideration,  and  proposes  a  R,C  damping 
circuit  added  to  the  system  to  obtain  stability. 

In  addition,  the  analysis  also  shows  that  in  order  to  maintain  the  stability  of  the 
system,  the  current  gain  of  the  current  controller  must  be  smaller  than  one,  i.e.,  k(s)  = 
I0/Iin  <1.0.    In  the  analysis,  first,  we  use  an  ideal  current-controller,  i.e.,  k=real 

constant.  Then,  based  on  PSpice  simulation,  we  use  a  simple  current-controller  model 
which  is  valid  up  to  1.1  kHz  for  20  kHz  switching  frequency.  From  these  two  kinds  of 
analyses,  the  effect  of  the  bandwidth  limitation  of  the  inverter  on  the  system  stability  is 
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revealed.  If  the  switching  frequency  is  not  high  enough,  more  detailed  model  needs  to 
be  considered  [70]. 

4.2  Stability  Analysis  of  the  Regenerative  Load  System  with  Source  Inductance 
When  a  regenerative  load  is  connected  to  the  utility  source,  the  system  may 
oscillate  because  of  the  presence  of  the  utility  source  impedance,  Ls.  The  schematic  of 
the  regenerative  load  system  with  source  inductance,  Ls,  is  shown  in  Fig.  4.1.  The 
PSpice  simulation  results  are  shown  in  Fig.  4.2. 

The  occurrence  of  such  instability  is  the  motivation  for  the  investigation 
reported  here.  There  are  many  related  research  works  on  the  stability  of  distributed 
power  supply  system  [57-68].  The  nature  of  the  problem  can  be  understood  by 
consideration  of  the  equivalent  circuit  of  the  feedback  port  of  the  regenerative  load 
system  as  shown  in  Fig.  4.3.  For  concreteness,  we  consider  the  case  for  the  DUT 
without  loading.  Hence,  the  input  impedance  of  the  DUT  can  be  represented  by  an 
inductance  Zli  =  sLl,  where  Ll  is  the  magnetization  inductance  of  the  input 
transformer.  This  equivalent  circuit  can  be  further  reduced  to  that  shown  in  Fig.  4.4, 
where  k  is  the  gain  of  the  current  controller.  The  terminal  voltage  Vi(s),  i.e.,  V(7,3), 
can  be  obtained  as 
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Fig.  4.1  The  schematic  circuit  diagram  of  the  "regenerative  load  system"  with  utility 
source  inductance  Ls. 
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Fig.  4.2  Simulated  current  and  voltage  waveforms  for  the  regenerative  load  system  with 
utility  source  inductance  Ls  =0.02mh. 

(a)  Utility  source  voltage  vs  =  v(7,3) ;  utility  source  current  is  =  i(vms) 

(b)  DUT  input  current  iL  =  i(vm) 

(c)  Current  controller  current  icc  =  i(vmc) 

(d)  DUT  output  voltage  v0  =  v(l) ;  DUT  output  current  i0  =  i(vmle) 

(e)  Capacitor  (Cdc)  voltage  Vdc  =  v(71,FGnd) 
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Fig.  4.3  The  equivalent  circuit  of  the  feedback  port  of  the 
regenerative  load  system. 
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Fig.  4.4  The  reduced  equivalent  circuit  of  Figure  4.3. 


where 


T  = 


ZL|| 


"  SCi 
1-k 


1-k 


k  =  gain  of  the  current  controller  (inverter) 


As  discussed  in  references  57  and  66,  the  impedance  ratio  at  the  interface,  defined  as  T 
=  Zs/Zeq',  can  be  considered  as  the  loop-gain  of  the  system  and  can  be  used  to 
determine  the  stability. 

The  root-locus  of  the  characteristic  equation 


is  shown  in  Fig.  4.5  for  k  =  0  to  k  =  0.99,  where  Ls  -  2mH,  L,i  =  600mH,  rti  =  0.00 1Q, 
CI  =  lOuF.  Figure  4.5  clearly  shows  one  pair  of  roots  on  the  imaginary  axis.  Hence, 
the  system  is  marginally  stable  and  has  sustained  oscillations.  The  simulated  results, 
using  Matlab,  is  shown  in  Fig.  4.6  for  k  =  0.99.  As  shown  in  Fig.  4.6,  the  oscillation 
frequency  is  opproximately  fos  =  U  KHZ.  This  is  in  agreement  with  f  =  a>/2n  = 
1 1 .5kHz,  which  is  calculated  from  the  root-locus  in  Fig.  4.5  for  k  =  0.99. 

If  the  gain  of  the  current  controller  is  increased  beyond  k  >1  where  the  root- 
locus  is  shown  in  Fig.  4.5(b),  we  can  then  see  that  there  is  one  pole  in  the  RHP.  The 
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Fig.  4.5  Root-locus  of  the  regenerative  load  system  with  source  impedance  Ls,  where 
Ls  =  2mH,  Li  =  600mH,  Ci=  lOuF.  (a)  stable  system,  k=0  to  0.99,+  for  k=0,  x  for 
k=0.99;  (b)  unstable  system,  k=1.001  to  1.10,  *  for  k=1.001,  o  fork=1.10. 
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Fig.  4.6  Matlab  simulation  results  of  the  regenerative  load  system  with  utility  source 
impedance  Ls,  where  Ls=2mH,  L/  =  600mH,  r/  =  0.00 1Q,  and  Ci  =10uF.  (a)  current 
of  capacitor  Ci;  (b)  voltage  at  the  interface  of  utility  source  and  the  DUT.  The 
oscillation  frequency  is  approximately,  fos  *  1.1  kHz  for  the  current  gain  k=0.99. 
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system  will  become  unstable  because  for  k  >1,  there  is  a  positive  feedback  with  loop- 
gain  greater  than  one. 

In  order  to  stabilize  the  system,  the  root-locus  must  be  shifted  from  the 
imaginary  axis  into  the  LHP.  Employment  of  state  feedback  and  pole  placement  can 
achieve  this  goal.  However,  the  simple  method  is  to  add  an  R,C  damping  circuit 
across  the  terminal  as  shown  in  Fig.  4.7.  The  Rp,Cp  damping  circuit  is  used  to  damp 
out  the  oscillation  for  the  case  of  current  gain  k  <  1 .  Therefore,  the  characteristic  roots 
of  the  system  will  move  into  the  LHP  as  shown  in  Fig.  4.8,  and  the  system  will 
become  asymtotically  stable  by  adding  this  damping  circuit  Rp,Cp. 

In  the  case  of  the  gain  of  the  current  controller  increasing  beyond  one,  k  >  1, 
there  will  be  one  pole  moved  into  the  RHP  as  shown  in  Fig.  4.9  and  the  system  will 
become  unstable  again.  The  reason  for  this  instability  is  due  to  the  positive  feedback 
with  loop-gain  greater  than  one,  k  >  1 .  For  the  stability  of  the  system,  the  damping 
circuit  must  be  added  and  at  the  same  time,  the  current  gain  of  the  current  controller 
need  to  be  smaller  than  one,  i.e.,  k  <  1 . 


4.3  Simulation  Results 
When  the  damping  circuit,  Rp,Cp,  is  added  to  the  terminal,  the  system  will 
become  asymptotically  stable  for  k  <  1,  as  the  root-locus  of  Fig.  4.8  (a)  shows.  The 
schematic  of  the  regenerative  load  system  with  the  damping  circuit  is  shown  in  Fig. 
4.10.    Figure  4.11  shows  the  simulation  results,  using  PSpice  software,  where 
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Fig.  4.7  The  equivalent  feedback  port  circuit  with  the  damping  circuit  included. 
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Fig.  4.8  Root-locus  of  the  regenerative  load  system  with  source  impedance  Ls,  and 
having  the  damping  circuit  Rp=10Q,Cp=10uF  added,  (a)  stable  system,  k=0  to 
k=0.99,  +  for  k=0,  x  for  k=0.99;  (b)  Zoom  of  (a)  around  the  origin  of  the  s  plane. 
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Fig.  4.9  Root-locus  of  the  regenerative  load  system  with  source  impedance  Ls,  and 
having  the  damping  circuit  Rp=10Q,  Cp=10uF  added,  (a)  unstable  system,  k=  1.001  to 
1 . 1 0,  *  for  k=  1 .00 1 ,  o  for  k=  1 . 1 0;  (b)  Zoom  of  (a)  near  the  imaginary  axis. 
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Fig.  4. 10  The  schematic  circuit  diagram  of  the  "regenerative  load  system"  with  the 
damping  citcuit  "R^Cp." 
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Fig.  4. 1 1  Simulated  current  and  voltage  waveforms  for  the  regenerative  load  system  with 
utility  source  inductance  Ls  =2mH,  after  the  damping  circuit  Rp=10  Q,  Cp=10uP,  is 
added. 

fa)  Utility  source  voltage  vs  =  v(7,3) ;  utility  source  current  is  =  i(vms) 

(b)  DUT  input  current  i  L  =  i(vm) 

(c)  Current  controller  current  icc  =  i(vmc) 

(d)  DUT  output  voltage  v0  =  v(l) ;  DUT  output  current  i0  =  i(vmle) 

(e)  Capacitor  (Cdc)  voltage  Vdc  =  v(71,FGnd) 

DUT  output  current  changes  from  a  peak  value  of  OA  to  10A  at  the  fifth  cycle  and  from 
1  OA  to  OA  at  the  tenth  cycle. 
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Ls=2mH,  Ci  =  lOuF  and  fs  =  25  kHz.  In  the  simulation  circuit,  the  DUT  output 
current  changes  from  peak  value  of  OA  to  10A  at  the  fifth  cycle  and  from  peak  value  of 
10A  to  OA  at  the  tenth  cycle. 


4.4  Discussion 

Up  to  this  point,  the  current  gain  of  the  inverter  is  assumed  as  a  real  constant. 
In  reality,  however,  the  current  gain  will  be  a  function  of  complex  frequency  s,  i.e., 
k(s),  as  derived  in  the  following.  The  equivalent  circuit  of  the  current  controller  is 
shown  in  Fig.  4.12.  From  Fig.  4.12,  the  output  current  IL2(s)  is  obtained  as 

.  kcl2[i;ef-kiscaieIL2]-V(7,3) 
IL2(s)  =   (4.3) 


where 


kcl2 


kc,  - 


kc2 


Therefore, 


sL2 

■ 

=  kcikc2 

Vdc  285 

Vtri  3 

:50 

e=ki=0.l(v/a) 

2.5mh 

v0 

L2  '    sL2  +kikcl2    sL2  +kjkcl2 


(4.4) 


If  the  term  (sL2  +  k,  k|2  )  is  very  large,  i.e.,  k  ki2  =  loop-gain  »  1,  then  the  second 
term  Vo  /  (sL2  +  k,  ki2  )  can  be  neglected,  and  thus 
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Iu(s)  sL2 
o — A  rcm-J  CONV. 


V(7,3) 


-    PWM  - 


4  O 


(l{(s)  +  Icl(s))-  KSaHa]a 


Fig.  4. 12  The  equivalent  circuit  of  the  unipolar  PWM  current  controller. 
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kjkci2 

T  ,  1 

(4.5) 


ML2-        L2         _  1 

1  ref        S  +  1  +  

L2  Wp 

where 

P  T 


For  I*ef  =  k0  [  II  (s)  +  Ici  (s)  ]  at  s=0,  the  actual  current  gain  will  be  given  as 

k(s)  =  ^  sk0—  (4.6) 

IL(s)  +  Icl(s)      °l  +  J_ 

where 

k0  =  dc  current  gain  of  the  inverter 
By  applying  k  (s)  of  Eq.  (4.6)  to  Eq.  (4.1),  the  terminal  voltage  of  the  system  without 
damping  circuit  is  obtained  as 


v,(s)  =  vs— ! 


1  +  T(s) 
_  y  a,s  +  a0 


(4.7) 


S  b3s3  +  b2s2  +  b,s+  b0 


where 

ai  =  L|,     ao  =  kj  kci2  U/L  2    \    b3  =  Ls  L|  C|, 
b2  =  Ls  L,C,  (1-ko)  kikci2/L2, 

bi  =  Ls  +  Li,      b0  =  k,  kci2  L|/L  2  +  (1-ko)  (kj  ken/  L  2)  Ls 


Hence,  the  characteristic  equation  of  the  system  without  damping  circuit  is 


b3S3  +  b2S2  +  b,S+bo  = 


(4.8) 


From  Eq.  (4.8),  the  root-locus  can  be  plotted  as  shown  in  Fig.  4.13.  In  Fig.  4.13,  there 
is  one  complex  pair  of  poles  in  the  RHP.  As  a  result,  the  system  will  be  unstable,  and 
the  oscillation  will  increase  with  time.  These  results  are  verified  in  Fig.  4.2  for  the 
actual  circuit  simulated  using  PSpice  software.  This  is  different  from  the  simulation 
results  in  Fig.  4.3,  which  assumes  k  =  real  constant,  and  the  oscillation  amplitude  will 
not  grow  with  time. 

Similarly,  the  system  with  the  damping  circuit  can  also  be  re-examined, 
assuming  the  current  gain  K(s)  of  Eq.  (4.6).  The  terminal  voltage  Vi(t)  will  be 
obtained  as 


V,(s)  =  Vs 


1  +  T(s) 


(4.9) 


1  + 


where 


damp      ^  p 


Therefore 
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Fig.  4.13  Root  locus  of  the  regenerative  load  system  with  utility  source  impedance 
Ls=2mH,  and  without  damping  circuit,  assuming  that  the  current  gain,  i.e., 
k(s)  =  k0  (kjkl2  /  L2)/(s+  kjk^  /  L2),  is  a  function  of  frequency,  where  k0  =  0  to 
1.02  ;  "+"  for  ko=0,  "x"  for  k0 =1.0  and  "o"  for  ko=1.02.  (a)  entire  root  locus;  (b) 
Zoom  of  (a)  near  the  imaginary  axis. 
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2 

Vfs)  =  V   a2S  +a's  +  ao   (4.10) 

1  sb4s4+b3s3+b2s2+b,s+b0 


where 

a2  =  RpCpL,,        ai  =  L,  (1  +  RpCp  kj  kci2/L2),      ao  =  Li  k,  kCi2/L2; 
b4  =  LsL,RpCpC,,    b3  -  LsL,  [  C,  +  Cp  +  (l-k0)RpCp(k,  kci2/L  2)  C,  ], 
b2  =  Ls  {L,  [  (l-kO)(  ki  kc,2/L2)  C,  +  (kj  kCi2/L2)  Cp]  +  RpCp}  +  RpCpL,, 
b,  =  Ls  {  1  +  (1-kO)  RpCp(ki  kc12/L2) }  +  L,  [  1  +  RpCp(k,  kci2/L2)  ], 
b0  -  (1-kO)  Ls  (kj  kci2/L2)  +  U  (ki  kci2/L2) 

and  the  characteristic  equation  of  the  system  with  damping  circuit  is  given  as 

b4S4  +  b3S3  +  b2S2  +  b,S  +  bo  =  0  (4-11) 

The  root-locus  of  the  system  can  be  obtained  from  Eq.  (4.11)  as  shown  in  Fig.  4.14. 
When  the  root-locus  in  Fig.  4.14  is  compared  with  that  of  Fig.  4.8  and  Fig.  4.9,  which 
assume  current  gain  k  =  real  constant,  it  can  be  seen  that  the  complex  dominant  pole 
pairs  in  Fig.  4.8  (b)  and  Fig.  4.14  (b)  are  almost  the  same,  and,  thus,  are  related  to  the 
damping  circuit  Rp,Cp.  The  difference  between  these  two  root-locus  is  in  the  other 
pole  pairs,  which  are  influenced  by  the  current  gain  k(s)  of  the  inverter. 

We  can  see,  however,  from  Fig.  4.14  (a),  that  the  system  is  stable  for  current 
gain  k0<  1 .  Therefore,  the  requirement  for  the  system  stability  is  that  the  dc  current 


74 


Root  Locus 


(a) 


Root  Locus 


(b) 


8000r 


imag 
-200fl- 


Fig.  4.14  Root  locus  of  the  regenerative  load  system  with  utility  source  impedance 
Ls=2mH,  and  with  damping  circuit,  Rp  =  10Q ,  Cp  =  lOuP ,  assuming  that  the  current 

gain,  i.e.,  k(s)  =  k0  (kjk12  /  L2)/(s+ kjk12  /L2),  is  a  function  of  frequency,  where 
k0  =  0  to  1.02;  "+"  for  ko=0,  "x"  for  ko=1.0  and  "o"  for  ko=1.02.  (a)  entire  root 
locus;  (b)  Zoom  of  (a)  near  the  imaginary  axis. 
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gain  k0  of  the  inverter  is  less  than  one,  i.e.,  k0  <  1,  and  the  damping  circuit  Rp,Cp  is 
added  to  the  regenerative  load  system. 

4.5  Summary 

This  chapter  presented  an  analysis  of  the  stability  of  the  regenerative  load 
system  with  the  utility  source  inductance,  Ls,  taken  into  consideration.  It  was  shown 
that  the  system  will  be  exponentially  unstable,  using  PSpice  simulation,  i.e.,  assuming 
that  the  current  gain  k(s)  of  the  current  controller  is  a  function  of  complex  frequency  s. 
If  the  current  gain  is  assumed  as  a  constant,  for  k  <  1,  the  system  is  still  marginally 
stable  with  sustained  oscillation,  as  shown  in  the  results  of  simulation  using  Matlab 
software.  These  instability  phenomena  are  corroborated  by  the  root-locus  analysis  of 
the  regenerative  load  system. 

Although  employment  of  state  feedback  and  pole  placement  can  relocate  these 
RHP  poles  into  the  LHP,  the  simple  method  is  to  add  an  R,C  damping  circuit  across 
the  interconnection  terminal  of  the  utility  source  and  the  DUT.  This  R,C  circuit  can 
damp  out  the  oscillation  for  the  case  of  current  gain  k  <  1.  Therefore,  the 
characteristic  roots  of  the  system  will  move  into  the  LHP  as  shown  in  the  root-locus 
analysis.  The  effect  of  the  R,C  damping  circuit  on  the  stability  of  the  system  is 
verified  using  PSpice  simulation. 

Even  with  the  addition  of  this  damping  circuit,  the  system  will  still  become 
unstable  if  the  current  gain  of  the  current  controller  is  greater  than  one,  i.e.,  k  >  1 . 


This  instability  is  demonstrated  by  the  movement  of  poles  into  the  RHP,  using 

the  root-locus  analysis.  The  physical  explanation  for  the  instability  is  that  if  the 
current  gain,  k,  of  the  current  controller  is  greater  than  one,  then  it  is  the  situation  of 
positive  feedback  with  loop-gain  greater  than  one  that  will  cause  instability.  In 
summary,  for  the  stability  of  the  system,  the  damping  circuit  must  be  added  across  the 
utility  source  terminal,  and  at  the  same  time,  the  current  gain  of  the  current  controller 
must  be  kept  below  one,  i.e.,|k(s)|  <  1 . 


CHAPTER  5 

INPUT  IMPEDANCE  SPECIFICATION 
OF  DEVICE  UNDER  TEST  FOR 
STABLE  REGENERATIVE  LOAD  SYSTEM 

5.1  Introduction 

In  order  to  find  the  input  impedance  specification  of  DUT  for  the  stability  of 
the  system,  there  are  many  related  methods  in  the  literature  [57-68].  Among  these 
methods,  the  two  subsystems  interaction  theorem  is  an  adequate  one  and  is  adopted 
here  [57,66].  First,  the  regenerative  load  system  is  transformed  into  two  equivalent 
subsystems;  then  the  two  subsystems  interaction  theorem  can  be  applied  to  analyze  the 
stability  of  the  system.  This  provides  an  alternative  analysis  method  in  addition  to  the 
root-locus  used  in  Chapter  4.  In  addition,  the  input  impedance  specification  of  the 
DUT  is  also  presented,  using  the  subsystem  interaction  theorem  [57,66]. 

At  first,  in  Sec.  5.2,  the  equivalent  impedance  concept  is  proposed  to  transform 
the  impedance  and  its  associated  controlled  current  source  into  an  equivalent 

Z 

impedance  such  that  Zeq  =  where  Z  is  the  original  impedance  and  k  is  the  gain 

(1  —  k) 

of  the  current  controller.  After  this  transformation,  the  system  can  be  further  divided 
into  two  equivalent  subsystems.  Then,  in  Sec.  5.4,  the  two  subsystems  interaction 
theorem  is  applied  to  analyze  the  stability  of  the  system.  We  found  that  the  system  has 
potential  stability  problems  for  the  active  DUT  load  (e.g.,  regulated  power  converters) 


77 


78 

due  to  the  phase  characteristics  of  the  small  signal  dynamic  input  impedance  of  the 
active  DUT  (negative  resistance  in  the  low-to-medium  frequency  range).  In  order  to 
clarify  the  effect  of  bandwidth  limitation  of  the  current  controller  on  the  system 
stability,  we  analyze  the  system  using  both  the  ideal  model,  i.e.,  current  gain  k=real 
constant,  and  the  simple  frequency-dependent  model  which  is  valid  up  to  1.1  kHz  for 
20  kHz  switching  frequency,  based  on  PSpice  simulation.  If  the  switching  frequency 
is  not  high  enough,  more  detailed  model  by  Tan  and  Middlebrook  [70]  needs  to  be 
considered.  From  the  stability  analysis,  we  found  that  the  bandwidth  limitation  of  the 
current  controller  has  an  adverse  effect  on  the  system  stability.  Therefore,  a  two 
damping  circuits  concept  is  proposed  in  Sec.  5.5  to  improve  the  stability  margin  of  the 
system.  With  the  addition  of  a  high  frequency  damping  circuit,  the  stability  margin 
increases  remarkably.  Subsequently,  in  Sec.  5.6,  the  input  impedance  specification  of 
the  DUT  for  a  specified  stability  margin  is  presented.  From  this  input  impedance 
specification  of  DUT,  we  concluded  in  Sec.  5.7  that  the  proposed  regenerative  load 
system  will  be  stable  for  the  practical  DUT  and  even  for  the  active  DUT.  The  only 
requirement  of  the  input  impedance  of  DUT  is  that  its  magnitude  is  greater  than  0.1Q 
for  all  the  frequency  range  of  interest.  The  validity  of  this  study  is  verified  by  PSpice 
simulation. 
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5.2  Defining  the  Regenerative  Load  System 
as  Two  Cascaded  Subsystems 


In  order  to  find  the  input  impedance  specification  of  the  DUT  for  a  stable 
regenerative  load  system,  the  utility  source  port  can  be  transformed  into  two  equivalent 
subsystems  as  shown  in  Fig.  5.1. 

Figure  5.1(a)  shows  the  original  input  port  equivalent  circuit  with  the  inverter 
represented  as  a  controlled  current  source.  This  controlled  current  source  can  be  split 
into  two  controlled  current  sources  as  shown  in  Fig.  5.1(b).  One  controlled  current 
source  is  parallel  with  the  capacitor  Ci,  and  the  other  is  parallel  with  the  input 
impedance  of  the  DUT. 

Since  the  controlled  current  source,  k  I/,  supplies  current  to  the  impedance  Z,, 
the  input  current  to  this  parallel  combination  will  be 


I/t>£?=I,-kI,  =(l-k)I, 


(5.1.) 


Hence  the  equivalent  load  impedance  is  given  as 


Z 


V 


v, 


l.eq  ~ 


I 


(l-k)I, 


l,eq 


z 


(5.2) 


1-k 


where 


input  impedance  of  DUT 


k=  current  gain  of  the  controlled  current  source  (inverter) 
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Fig.  5.1  The  transformed  equivalent  subsystems  of  the  regenerative  load  system. 
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Fig.  5.1  Continued. 


Similarly,  the  parallel  combination  of  the  l/(sCi)  and  klci  can  be  transformed 
into  an  equivalent  impedance,  Zci,eq,  as 

Zc.eq=^^  (5-3) 

This  transformed  equivalent  system  can  be  further  divided  at  the  terminal  A,  B  as 
shown  in  Fig.  5.1(c).  The  left  part  will  be  represented  by  a  Thevenin's  equivalent 
circuit  with  the  Vs,eq  and  Zs,eq  given  as 

Vs,eq  =  Vs     Zc'eq  (5.4) 

+  ^cl,eq 

IZjZqJ  gq 
^cl,eq  =  ~Z      j  (5-5) 

When  the  DUT  is  connected  to  the  Thevenin's  equivalent  source  circuit,  the 
regenerative  load  system  will  be  reduced  to  two  cascaded  subsystems  as  shown  in  Fig. 
5.1(d).  These  two  subsystems  are  the  (1)  equivalent  source  subsystem  and  (2) 
equivalent  DUT  load  subsystem. 

In  the  next  section,  a  methodology  for  developing  the  input  impedance 
specification  of  an  unknown  DUT  subsystem  for  stable  system  operation  is  introduced. 
The  method  assumes  that  the  output  impedance  of  the  utility  source  subsystem  is 
known.  If  the  DUT  meets  the  input  impedance  specification,  the  stability  of  the 
regenerative  load  system  will  be  ensured,  despite  allowing  impedance  overlap  at  the 
interface  of  the  subsystems  [57,66]. 
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5.3  System  Interaction  Analysis 
5.3.1  Impedance  Ratio  as  a  Loop-Gain 

Consider  the  integration  of  the  two  subsystems  as  shown  in  Fig.  5.2. 
Assuming  that  each  subsystem  is  individually  stable,  the  overall  transfer  function  of 
the  integrated  system  can  be  expressed  as  the  product  of  the  individual  subsystem 
transfer  functions  (  Fa  and  Fb  )  and  a  loading  factor  [57-68]. 


where 


V2B  _  FAFB 

7. 


FaFb 


1  +  Tm 


(5.6) 


Z 

Tm  =      =  Equivalent  loop-gain  of  the  integrated  system 


Zin 


FA  =  Forward  transfer  function  of  load  subsystem 

without  load  impedance  ,  i.e.,  Zl=  <» 
FB  =  Forward  transfer  function  of  load  subsystem 

without  source  impedance ,  i.e.,Zs=  0 
Zo=  Output  impedance  of  the  source  subsystem 
Zin=  Input  impedance  of  the  load  subsystem 

1      „      Zin  . 

1  +  Tm  Z0+Zin 


Loading  factor 
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Fig.  5.2  Two  cascaded  subsystems. 


Both  Zo  and  FA  are  input-to-output  transfer  functions  of  the  source  subsystem 
and  have  common  denominators  (characteristic  equation  of  the  source  subsystem). 
Also,  Yin  (  1/Zin  )  and  FB  are  input-to-output  transfer  functions  of  the  load  subsystem 
and  have  common  denominators  (characteristic  equation  of  the  load  subsystem).  The 
above  statement  is  true  when  both  the  subsystems  are  observable  at  the  interface  [71]. 

Since  the  denominator  of  FAFB  is  canceled  by  the  numerator  of  the  loading 
factor  in  Eq.  (5.6)  (also  the  denominator  of  Tm),  the  characteristic  equation  of  the 
integrated  system  will  be 

1  +  Tm=0  (5.7) 

As  discussed  in  references  58  and  66,  the  impedance  ratio  at  the  interface,  defined  as 
Tm=Zo/Zin,  can  be  viewed  as  the  equivalent  loop-gain  of  the  system.  This  equivalent 
loop-gain  can  be  used  to  determine  the  stability  and  performance  of  the  system. 

The  stability  of  the  integrated  system  can  be  determined  by  calculating  the 
zeros  of  Eq.  (5.7),  i.e.,  root-locus  plot,  or  by  applying  the  Nyquist  criteria  to  this  loop- 
gain  Tm.  The  root-locus  method  is  adopted  in  Chapter  4  and  in  this  chapter,  the 
Nyquist  criteria  is  used.  The  root-locus  method  requires  analytical  models;  however, 
the  Nyquist  plot  can  easily  be  determined  from  either  analytical  models  or 
measurement  data.  Also,  it  should  be  noted  that  the  choice  of  source  and  load 
subsystems  is  not  dependent  on  the  direction  of  power  flow.  However,  the  source 
subsystem  must  be  chosen  such  that  its  output  impedance  at  the  interface  bus  is  an 
input-to-output  transfer  function  of  the  source  subsystem.     Similarly,  the  load 


subsystem  must  be  chosen  such  that  its  input  admittance  at  the  interface  bus  is  an 
input-to-output  transfer  function  of  the  load  subsystem. 

If  |Zin|»|Z0|  for  all  frequencies,  the  loading  effect  is  negligible  and  the 
system  stability  will  depend  only  upon  the  stability  of  the  individual  subsystems. 
When  |Z0|  is  larger  than  |Zin|,  a  considerable  loading  effect  exists.    However,  a 

loading  effect  does  not  necessarily  imply  a  stability  problem.  In  this  case,  the  Nyquist 
criteria  can  be  applied  to  the  loop-gain,  Tm,  to  determine  the  system  stability. 

5.3.2  Effect  of  Impedance  Overlap  on  System  Stability  and  Performance 

For  the  case  where  |Z0|  exceeds  |Zin|  for  some  frequencies,  further  analysis  is 

needed  to  determine  the  system  stability  as  discussed  in  references  57  and  66.  An 
example  of  two  impedance  overlapping  is  shown  in  Fig.  5.3.  In  the  frequency  range 
where  I  Zin  I  <  I  Zo  I ,  the  magnitude  of  loop-gain  Tm  is  greater  than  0  dB,  as  shown  in 
Fig.  5.3(a),  (b),  (c).  The  magnitude  of  Tm  crosses  the  0  dB  line  at  the  two  frequencies 
where  I  Zin  I  =  I  Zo  I  [points  A  and  B  in  Fig.  5.3  (a),  (b),  (c)].  In  addition,  the 
magnitude  of  Tm  must  approach  zero  as  the  angular  frequency,  (D=2nf,  becomes 
infinity  for  the  applicability  of  the  Nyquist  criteria. 


87 


Fig.  5.3  Impedance  overlap  at  the  interface,  (a)  impedance  comparison;  (b)  phase  plot 
of  impedance;  (c)  Bode  magnitude  plot  of  Tm  =  Z0/Zin  ;  (d)  Bode  phase  plot  of  Tm  ; 
(e)  polar  plot  of  Tm  =  Z0/Zin  (on  the  next  page). 


Fig.  5.3  Continued,  (e)  polar  plot  of  Tm  =  Z0/Zj 
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Since  the  |Tm|  crosses  the  0  dB  line  twice,  two  different  phase  margins  can  be 

defined  assuming  the  subsystems  are  individually  stable,  one  at  the  initial  point  of  the 
overlap  [point  A  in  Fig.  5.3  (a),  (b),  and  (c)]  and  the  other  at  the  final  point  of  the 
overlap  [point  B  in  Fig.  5.3(a),  (b),  (c)].  For  this  example,  the  phase  margin  at  the 
initial  point  of  overlap,  PMi  ,  is  given  by 

PM,=180°-ZTm  (5.8) 

and  the  phase  margin  at  the  final  point  of  overlap,  PM2 ,  is  given  by 

PM2=180°+ZTm  (5.9) 

Figure  5.3(d)  shows  the  polar  plot  of  Tm  together  with  two  phase  margins. 

In  addition  to  the  stability  of  the  loop-gain,  Tm,  the  phase  margins  of  Tm  are 
also  of  particular  importance  in  the  cascading  of  two  subsystems.  A  small  phase 
margin  drastically  increases  the  magnitude  of  the  loading  factor  1/(1+Tm),  resulting  in 
severe  peaking  in  the  closed-loop  transfer  function  of  the  integrated  system.  As  a 

conservative  design  practice,  a  phase  margin  of  60°  can  be  applied  for  most  practical 
applications. 
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5.4  Stability  Analysis  of  the  Regenerative  Load  System 
as  Two  Cascaded  Subsystems 

5.4. 1 .  Stability  of  Subsystems 

Section  5.3  mentioned  that  if  the  two  subsystems  are  individually  stable,  the 
impedance  ratio  at  the  interface  as  in  Eq.  (5.6)  can  be  considered  as  the  loop-gain  of 
the  integrated  system,  and  the  Nyquist  criterion  can  be  applied  to  the  loop-gain,  Tm,  as 
defined  in  Eq.  (5.6)  to  determine  the  system  stability.  Therefore,  the  source  subsystem 
needs  to  be  stabilized  by  adding  a  damping  circuit  Rp,Cp  as  shown  in  Fig.  5.4.  In 
addition,  this  compensated  regenerative  load  system  can  be  transformed  into  two 
equivalent  subsystems,  as  shown  in  Fig.  5.4(c),  (d),  and  (e)  where 


VTH  -  Vs 


ZP 

Zcl.eq 

zs  +  (zp 

Zcl.eq  ) 

(5.10) 


Z-th  -  zs|zp 


'cl.eq 


(5.H) 


and  where 


Jcl.eq 


J  1_ 

sC,  1-k 


z>=R>+A 


Zs  =  sLs 


The  stability  of  the  compensated  source  subsystem,  where  Rp=10Q  and 
Cp=10uF,  is  investigated  by  the  root-locus  plot  of  the  source  subsystem,  using  the 
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Fig.  5.4  The  transformed  equivalent  subsystems  of  the  compensated  regenerative  load 
system. 


(c)  two  subsystems  with 
equivalent  impedance: 


(d)  two  cascaded  subsystems 


o-  -o- 


(e)  represented  as  Thevenin's  equivalent  circuit. 


Fig.  5.4  Continued. 
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current  gain  K(s)  as  the  varying  parameter,  where  the  current  gain  is  given  as  in  Eq. 


(4.6) 


and  where 


K(s)  =  k0 


kjki2 
s+k|kl2 


L2  (5.12) 
1 


1  + 


p 


k0  =  dc  current  gain  of  the  controlled  current  source  (inverter) 
kj  =0.1  [V/A] 


k12-  kjk2 


V 

vtn 


k2=50 

L2  =2.5mH=  inductance  of  the  current  controller 
kjkj2 

con  = 


The  root-locus  in  Fig.  5.5  shows  that  if  the  dc  current  gain,  i.e.,  k0  in  Eq.  5.12  is 

greater  than  one,  then  the  root-locus  will  enter  into  the  RHP  (indicated  as  point  o),  and 
the  source  subsystem  will  become  unstable.  Hence,  the  dc  current  gain,  k0,  of  the 

inverter  must  remain  smaller  than  one,  i.e.,  k0<  1,  for  stability.   Figure  5.6  shows 
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Fig.  5.5  Root-locus  of  the  source  subsystem  as  the  function  of  current 
gain,k(s)  =  k0  (kjk12  /  L2)/(s+  kjk|2  /  L2),  where  k0  from  0  to  1.10,  and  damping 
circuit  Rp=10Q,  Cp=10uP.  (a)  entire  root  locus;  (b)  zoom  of  (a)  near  the  imaginary 
axis;  (c)  Zoom  of  (b)  around  the  origin  (on  the  next  page). 
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Fig.  5.5  Continued,  (c)  zoom  of  (b)  around  the  origin. 
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Fig.  5.6  Root  locus  of  source  subsystem  as  function  of  damping  resistor  Rp,where  the 
current  gain  k0  =  0.95.  (a)  entire  root  locus;  (b)  Zoom  of  (a)  near  the  imaginary  axis; 
(c)  zoom  of  (b)  around  the  origin  (on  the  next  page). 


Fig.  5.6  Continued,  (c)  zoom  of  (b)  around  the  origin. 
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another  root-locus  of  the  source  subsystem,  where  dc  current  gain  k0-0.95,  and  the 
damping  resistor  Rp  is  used  as  a  varying  parameter.  Figure  5.6  shows  there  are  a  pair 
of  poles  on  the  imaginary  axis  when  Rp=0  as  shown  in  Fig.  5.6(c),  marked  as  points  +. 
However,  as  Rp  increases  from  zero  to  Rp=50Q,  the  other  pair  of  poles  will  move 
from  the  far  end  of  LHP  towards  the  imaginary  axis  as  shown  in  Fig.  5.6(a)  and  (b). 
Based  on  these  observations,  the  value  of  the  resistance,  Rp=10Q,  is  chosen  as  the 
adequate  one.  This  choice  has  been  justified  by  the  PSpice  simulation  of  the  source 
subsystem. 

The  DUT  subsystem  in  Fig.  5.4  can  be  divided  into  two  modes,  nonconduction 
and  conduction,  as  shown  in  Fig.  5.7(a)  and  (b),  depending  on  the  conduction 
conditions  of  the  bridge  rectifier  of  Fig.  3.3  in  Chapter  3.  In  Fig.  5.7(a),  the  L/i 
represents  the  leakage  inductance,  and  the  L/t  is  the  magnetization  inductance  of  the 
transformer.  The  r/i  and  Vr  are  the  coil  resistance  and  equivalent  resistance  of  core 
losses,  respectively.  The  C  and  Rl  of  Fig.  5.7(b)  represent  the  passive  DUT  load  in 
conduction  mode.  Figure  5.7(c)  shows  the  equivalent  small  signal  DUT  subsystem, 
where  the  negative  resistance,  Rn=  -10Q,  represents  the  worst  case  small  signal 
impedance  of  the  dc-to-dc  regulator  [9].  The  Rb,Cb  circuit  in  Fig.  5.7(c)  is  used  to 
damp  the  oscillation  caused  by  the  negative  dynamic  resistance  Rn  of  the  regulated  dc- 
to-dc  converter.  The  stability  of  all  the  modes  in  Fig.  5.7  have  been  verified  through 
the  PSpice  simulation  where  k=0.99  (or  k0=0.99). 
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Fig.  5.7  The  equivalent  DUT  load  subsystems. 
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5.4.2  Stability  of  the  Integrated  System 

The  stability  of  the  regenerative  load  system  can  be  investigated  by  cascading 
the  three  modes  of  DUT  subsystems  of  Fig.  5.7  to  the  source  subsystem  as  shown  in 
Fig.  5.4.  The  stability  analysis  for  the  non-conduction  mode  of  DUT  is  shown  in  Fig. 
5.8  where  r/i=0.001Q  and  L/i=25uH  are  the  internal  resistance  and  leakage  inductance 
of  the  transformer,  and  L/t=600mH  and  r/t=0.01Q  are  the  magnetization  inductance  and 
equivalent  resistance  of  core  losses.  Figure  5.7(a)  represents  an  equivalent  circuit  of 
the  transformer  with  its  secondary  open  circuited. 

In  Fig.  5.8(a),  the  current  gain  is  assumed  as  a  real  constant,  i.e.,  k=0.99.  We 
can  see  I  Zo  I «  I  Zin  I  for  all  frequencies,  and  therefore,  there  is  no  impedance 
overlap  and  the  loop-gain  |  Tm  | «  1  from  co  «  0  to  co  «  qo.    Since  there  is  no 
encirclement  of  the  (-1,0)  point  in  the  polar  plot  of  Tm  in  Fig.  5.8(d),  the  integrated 

system  will  be  stable  by  applying  the  Nyquist  criterion. 

Fig.  5.9  shows  the  same  analysis  as  in  Fig.  5.8;  the  current  gain,  however,  is 
assumed  as  a  function  of  frequency,  i.e.,  K(s)  as  in  Eq.  (5.12).  Figures  5.9(a)  and  (d), 
clearly  show  that  there  is  no  impedance  overlap  and  the  equivalent  loop-gain  |  Tm  | « 
1.  Hence  the  polar  plot  of  Tm  will  not  encircle  the  (-1,0)  point,  and  the  system  will 
be  stable. 

For  the  conduction  mode  with  passive  load  resistance  as  in  Fig.  5.7(b),  the 
results  of  stability  analysis  of  the  regenerative  load  system  are  shown  in  Fig.  5.10, 
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Fig.  5.8  Stability  analysis  of  the  regenerative  load  system  for  the  nonconduction  mode 
of  DUT,  i.e.  Fig.  5.7(a),  and  the  current  gain  is  assumed  as  k=0.99(real  constant),  (a) 
impedance  comparison;  (b)  phase  plot  of  impedance;  (c)  Bode  magnitude  plot  of 
Tm  =  Z0  /Zm  ;  (d)  Bode  phase  plot  of  Tm ;  (e)  polar  plot  of  Tm  (on  the  next  page). 
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Fig.  5.8  Continued,  (e)  polar  plot  of  T, 
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Fig.  5.9  Stability  analysis  of  the  regenerative  load  system  in  the  non-conduction  mode 
DUT  of  Fig.  5.7(a),  assuming  k(s)  =  k0  (kjk12  /  L2)/(s  +  kjk12  /  L2),  and 
k0  =0.99.  (a)  impedance  comparison;  (b)  phase  plot  of  impedance;  (c)  Bode 
magnitude  plot  of  Tm  =  Z0  /  Zin ;  (d)  Bode  phase  plot  of  Tm ;  and  (e)  polar  plot  of  Tm 
(on  the  next  page). 


Fig.  5.9  Continued,  (e)  polar  plot  of  T, 
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Fig.  5.10  Stability  analysis  for  the  conduction  mode  of  Fig.  5.7(b),  where  the  current 
gain  k=0.99  (real  constant),  (a)  impedance  comparison;  (b)  phase  plot  of  impedance; 
(c)Bode  magnitude  plot  of  Tm  =  Z0/Zin  ;  (d)  Bode  phase  plot  of  Tm  ;  (e)  polar  plot  of 
Tm ;  and  (f)  zoom  of  (e)  around  the  origin  (on  the  next  page). 


Fig.  5.10  Continued,  (e)  polar  plot  of  Tm ;  and  (f)  zoom  of  (e)  around  the  origin. 
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where  current  gain  is  assumed  as  a  real  constant  of  value  k=0.99.  In  Fig.  5.10(a),  there 
is  impedance  overlap  between  points  A  and  B,  such  that  the  equivalent  loop-gain, 
Tm=Zo/Zin,  will  be  greater  than  one  in  the  frequency  interval  from  A  to  B.  However, 
there  are  phase  margins  PM|  and  PM2  as  can  be  seen  in  Fig.  5.10(d)  and  (e).  The  polar 
plot  will  not  encircle  the  (-1,0)  point  as  shown  in  Fig.  5.10(e)  and  therefore,  the 

regenerative  load  system  will  be  stable. 

In  Fig.  5.1 1,  the  same  conduction  mode  of  Fig.  5.7(b)  is  analyzed,  assuming  the 
current  gain  K(s)  as  in  Eq.  (5.12),  where  the  resistance  Rp  is  changed  to  2Q  in  order  to 
make  the  system  stable.  Fig.  5.1 1(a)  shows  two  times  of  impedance  overlaps.  Hence, 
the  equivalent  loop-gain  Tm  will,  be  greater  than  one  over  two  intervals  as  shown  in 
Fig.  5.1 1(c),  (e),  and  (f).  For  the  second  interval  of  impedance  overlap,  the  equivalent 
loop  gain  Tm  will  almost  contact  the  (-1,0)  point  as  shown  (point  C)  in  Fig.  5.1 1(f). 

Therefore,  the  regenerative  load  system  will  become  only  marginally  stable. 

If  the  DUT  is  a  regulated  dc-to-dc  converter,  the  equivalent  small  signal 
dynamic  input  impedance  will  be  a  negative  resistance  Rn  as  in  Fig.  5.7(c).  The  above 
negative  load  resistance  can  be  considered  as  a  worst-case  when  performing  subsystem 
interaction  analysis  as  done  by  Lewis  et  al.  [63].  The  results  of  stability  analysis  of 
two  subsystems  interaction  with  negative  dynamic  load  resistance  is  shown  in  Fig. 
5.12  where  the  current  gain  of  the  inverter  is  assumed  as  a  constant  of  k=0.99.  Figure 
5.12  shows  only  one  interval  of  impedance  overlap  between  points  A  and  B,  and  the 
polar  plot  of  equivalent  loop-gain,  Tm=Zo/Zin  as  in  Fig.  5.12(e),  shows  there  are  large 
phase  margins  PMi  and  PM2.  Therefore,  the  regenerative  load  system  will  be  stable. 
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Fig.  5.11  Stability  analysis  for  the  conduction  mode  of  Fig.  5.7(b),  assuming  current 
gain  is  a  function  of  frequency,  i.e.,  k(s)  =  k0  (kjk12  / L2)/(s+  kjk^  / L2)  and 
k0  =0.99.  Rp  is  2  Q.  (a)  impedance  comparison;  (b)  phase  plot  of  impedance; 
(c)Bode  magnitude  plot  of  Tm  =  Z0/Zin  ;  (d)  Bode  phase  plot  of  Tm ;  (e)  polar  plot  of 
Tm ;  and  (f)  zoom  of  (e)  around  the  origin  (on  the  next  page). 


Fig.  5.11  Continued,  (e)  polar  plot  of  Tm ;  and  (f)  zoom  of  (e)  around  the  origin. 
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However,  if  the  current  gain  is  assumed  as  frequency  dependence,  i.e.,  K(s),  as 
in  Eq.  (5.12),  there  will  be  two  impedance  overlaps  as  shown  in  Fig.  5.13(a)  where  the 
Rp  is  equal  to  2Q.  The  polar  plot  of  equivalent  loop-gain  as  in  Fig.  5.13(e)  will  have 
two  loops  intersected  with  the  unit  circle.  The  second  loop  intersects  the  unit  circle  at 
point  C  as  shown  in  Fig.  5.13(e)  with  a  small  phase  margin.  Therefore,  the 
regenerative  load  system  will  have  only  marginal  stability  if  current  gain  is  frequency 
dependent  as  in  Eq.  (5.13). 

Comparing  Fig.  5.13(a)  with  Fig.  5.12(a),  we  see  that  the  second  impedance 

overlap  interval  occurs  approximately  at  frequency  m  2  x  104  Hz.  From  Fig.  5.4  and 
Fig.  5.7,  we  can  see  that  both  the  Zo  and  Zin  are  functions  of  current  gain  k(s).  The 
Bode  plot  of  k(s)  is  shown  in  Fig.  5.14.  As  shown  in  Fig.  5.14(b),  the  phase  of  current 

gain  k(s)  has  appreciable  lag  around  frequency  *  2  x  104Hz.  This  phase  lag  causes 

I  Zin  I  to  decrease  and  I  Zo  I  to  peak  around  f  *  2  x  104Hz  as  shown  in  Fig.  5.13(a). 
This  impedance  overlap  causes  another  loop  of  Tm  with  magnitude  greater  than  one, 
as  shown  in  Fig.  5.13(e).  This  loop  intersects  the  unit  circle  at  point  C,  near  the  (-1,0) 

point,  thus  the  system  is  marginally  stable.  We  concluded  that  the  bandwidth 
limitation  or  non-ideality  of  current  gain  k(s)  at  high  frequency  will  deteriorate  the 
stability  of  the  integrated  system. 


Fig.  5.12  Stability  analysis  for  the  conduction  mode  of  Fig.  5.7(c),  i.e.,  negative  small 
signal  dynamic  resistance  and  assuming  current  gain  k=0.99  (real  constant),  (a) 
impedance  comparison;  (b)  phase  plot  of  the  impedance;  (c)  Bode  magnitude  plot  of 
Tm  =  Z0  /Zin  ;  (d)  bode  phase  plot  of  Tm  ;  (e)  polar  plot  of  Tm  (on  the  next  page). 


Fig.  5.12  Continued,  (e)  polar  plot  of  T, 
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Fig.  5.13  Stability  analysis  for  the  conduction  mode  of  Fig.  5.7(c),  i.e.,  negative  small 
signal  dynamic  resistance,  and  assuming  the  current  gain  is  a  function  of  freqency,  i.e. 
k(s)  =  k0 (kik12  /L2)/(s+kik12  /L2),  where  k0  =0.99(real  constant).  Rp  is  2  Q. 
(a)  impedance  comparison;  (b)  phase  plot  of  impedance;  (c)  magnitude  plot  of  Tm ;  (d) 
phase  plot  of  Tm ;  (e)  polar  plot  of  Tm  (on  the  next  page). 
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5.5  Stability  Improvement  of  Regenerative  Load  System 

In  Section  5.4,  there  is  an  additional  impedance  overlap  around  f  =  2  x  104  Hz 
as  shown  in  Fig.  5.11(a)  and  Fig.  5.13(a)  due  to  the  non-ideality  of  the  current  gain 
k(s)  as  shown  in  Fig.  5.14.  This  additional  impedance  overlap  will  cause  another 
intersection  of  the  polar  plot  of  Tm  with  the  unit  circle  at  point  C  as  shown  in  Fig. 
5.1 1(e)  and  Fig.  5.13(e),  respectively.  The  phase  margin  at  point  C  is  very  small  and 
will  deteriorate  the  stability  of  the  regenerative  load  system. 

The  remedy  for  this  stability  deterioration  is  to  add  another  high  frequency 
damping  circuit  Rpi,Cpi  to  the  source  subsystem  as  shown  in  Fig.  5.15  such  that 


f3dB  =  0     *       *2xl04  Hz  (5.13) 


ZP1 


f=f,3dB 


=  V2Rpl  (5.14) 


If  Rpl  is  chosen  as  1  Q,  then  Zpl  »  1.414Q  and  Cp  *  8uF.   Actually,  if 

Cpl=2uF  and  Rpl=lQ  is  used  as  the  high  frequency  damping  circuit  for  the  source 
subsystem,  an  adequate  stability  improvement  of  the  regenerative  load  system  can  be 
obtained. 

Fig.  5.16(a)  shows  only  one  impedance  overlap  between  points  A  and  B 
compared  to  the  two  overlaps  in  Fig.  5. 1 1(a).  The  elimination  of  the  second  overlap  is 

due  to  the  lowering  of  |Zo  |  around  the  frequency  f  »  2  x  104  Hz  by  the  damping 
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(b)  Two  subsystems 


Fig.  5  .15  The  equivalent  subsystems  of  the  regenerative  load  system  with 
two  damping  circuits  RppCp,  and  Rp,Cp. 
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Fig.  5.16  Stability  analysis  for  the  conduction  mode  of  Fig.  5.7(b),  i.e.  passive 
resistance  DUT,  and  assuming  the  current  gain  is  a  function  of  frequency,  i.e. 
k(s)  =  k0(kikI2  /L2)/(s  +  kjk12  /L2),  where  ko=0.99.  Rp  is  2  Q  and  the  high 
frequency  damping  circuit,  Rpl  =1Q,  Cpl=2uF,  is  added,  (a)  impedance  comparison; 
(b)  phase  plot  of  impedance;  (c)  Bode  magnitude  plot  of  Tm ;  (d)  Bode  phase  plot  of 
Tm  ;  (e)  polar  plot  of  Tm  (on  the  next  page). 
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circuit  Rpi  and  Cpi.  The  polar  plot  of  the  loop-gain  Tm=Zo/Zin  has  only  one  loop  as 
shown  in  Fig.  5.16(e)  with  quite  large  phase  margins  PMi  and  PM2.  Compared  to  the 
polar  plot  of  Tm  in  Fig.  5.1 1(f),  Fig.  5.16(e)  shows  that  the  stability  of  the  regenerative 
load  system  has  been  improved  notably. 

Similarly,  Fig.  5.17  shows  the  stability  analysis  for  the  conduction  mode  of  the 
active  regulated  dc-to-dc  converter  with  the  high  frequency  damping  circuit  Rpi=lQ, 
Cpi=2uF  added  to  the  source  subsystem.  In  comparison  with  Fig.  5.13(a),  we  can  see 
that  there  is  no  impedance  overlap  in  Fig.  5.17(a).  Hence,  the  regenerative  load 
system  will  be  stable  as  shown  in  the  polar  plot  of  Fig.  5.17(e). 

5.6  Input  Impedance  Specification  of  Device  Under  Test 
We  have  seen  in  Sec.  5.4  and  5.5  that  by  applying  the  Nyquist  criteria  to  a 
polar  plot  of  the  loop-gain,  Tm,  we  can  examine  both  the  stability  and  relative 
stability.  In  this  section,  we  develop  a  methodology  for  deriving  the  DUT  input 
impedance  specification,  assuming  that  the  output  impedance  of  the  source  subsystem 
is  known.  The  method  ensures  desired  stability  margins  despite  allowing  impedance 
overlap  at  the  interface  between  subsystems. 

5.6.1  Forbidden  Region 

In  considering  Fig.  5.18,  suppose  that  the  loop  gain,  Tm,  never  enters  the 
shaded  portion  of  the  polar  plot.    In  such  a  case,  assuming  that  the  individual 
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Fig.  5.17  Stability  analysis  for  the  conduction  mode  of  Fig.  5.7(c),  i.e.,  negative  small 
signal  dynamic  resistance,  and  assuming  the  current  gain  is  function  of  frequency,  i.e., 
k(s)  =  k0(kik,2  /L2)/(s+kik|2  /L2),  where  ko=0.99.  Rp  is  2  Q  and  the  high 
frequency  damping  circuit,  Rpl  =1Q,  Cpl=2uF,  is  added,  (a)  impedance  comparison; 
(b)  phase  plot  of  impedance;  (c)  Bode  magnitude  plot  of  Tm ;  (d)  Bode  phase  plot  of 
Tm ;  and  (e)  polar  plot  of  Tm  (on  the  next  page). 
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subsystems  are  stable  and  all  eigenvalues  are  observable,  it  is  not  possible  for  the  polar 
plot  of  Tm  to  encircle  the  (-1,0)  point  and  so  the  integrated  system  stability  is 
guaranteed.  Also,  anytime  Tm  crosses  the  unit  circle,  the  phase  of  Tm  will  always  be 

such  that  there  is  at  least  60°  of  the  phase  margin.  Furthermore,  anytime  Tm  crosses 
the  negative  portion  of  the  real  axis,  the  magnitude  of  Tm  will  be  less  than  0.5, 
ensuring  at  least  6  dB  of  gain  margin. 

If  the  polar  plot  of  Tm  avoids  this  forbidden  region,  this  condition  guarantees 
that  the  integrated  system  will  be  stable,  will  not  be  conditionally  stable,  will  have  a 
desired  stability  margin,  and  will  have  minimal  performance  degradation. 

5.6.2  Formation  of  the  Load  Impedance  Specification 

By  using  the  output  impedance  curve  of  the  source  subsystem,  it  can  be  shown 
how  to  generate  an  impedance  specification  for  the  input  impedance  of  the  subsystem 
which  causes  Tm=Zo/Zin  to  avoid  the  forbidden  region  of  Fig.  5.18.  The  forbidden 
region  shown  in  Fig.  5.18  corresponds  to  the  shaded  area  determined  by  the  following 
two  inequalities 

|Z0|-|Zin|>-GM[dB]  (5.15) 
360°N  +  (180°  -PM,)<  ZZ0  -ZZin  <  360° N  +  (180°  +  PM2)  (5.16) 

where 

N=0,  ±1  ,  ±2,... 
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Based  on  Eq.  (5.15)  and  Eq.  (5.16),  two  important  concepts  of  gain  limit  and 
unacceptable  phase  band  can  be  developed. 


Gain  limit 

The  gain  limit,  |  Z  |  Limit ,  is  defined  as 

lZlLimit  =|Z0|  +  GM     [dB]  (5.17) 
When  |Zin|  stays  above  the  gain  limit  for  all  frequency,  the  polar  plot  does  not  enter 
the  forbidden  region,  and  the  load  subsystem  automatically  meets  the  input  impedance 
specification.  When  |Zin|  falls  below  the  gain  limit  for  some  frequency,  the  polar  plot 

of  Tm  could  enter  the  forbidden  region.  In  this  case,  the  phase  of  the  load  impedance 
should  be  tested. 


Unacceptable  phase  band 


The  unacceptable  phase  band  defines  the  region  that  the  ZZin  should  not  enter 
to  achieve  the  desired  phase  margin.  For  a  given  ZZ0  and  desired  phase  margin,  the 
unacceptable  phase  band  can  be  easily  obtained,  using  Eq.  (5.16)  as 

ZZ0  -[360°N  +  (180°  +PM2)]  <  ZZin>forbidden  < 

ZZo-[360°N  +  (180°-PM,)]  (5.18) 
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In  order  to  avoid  the  forbidden  region,  ZZin  should  not  enter  the  unacceptable  phase 
bands  for  all  frequencies  where  the  |Zin|  is  lower  than  the  gain  limit. 

5.6.3  Equivalent  Input  Impedance  Specification 

The  use  of  the  input  impedance  specification  applied  to  the  conduction  mode 

of  Fig.  5.7(b)  is  shown  in  Fig.  5.19,  assuming  GM=6dB  and  PM=50° .  Since  the  |Zin| 

is  lower  than  the  gain  limit  around  the  103  Hz  and  2  x  104  Hz  as  shown  in  Fig. 
5.19(a),  the  ZZin  must  be  checked  against  the  phase  specification.  As  shown  in  Fig. 
5.19(b),  the  phase  of  the  ZZin  enters  the  unacceptable  phase  bands  at  points  A  and  C, 
where  |Zin|  falls  below  the  magnitude  limit .  Therefore,  as  shown  in  Fig.  5.19(f),  the 
polar  plot  of  the  loop-gain  Tm  does  enter  the  forbidden  region,  although  the  integrated 
system  is  still  marginally  stable. 

If  the  high  frequency  damping  circuit  Rpl=lQ,  Cpl=2uF  is  added  to  the 
source  subsystem,  the  equivalent  DUT  input  impedance  specification  becomes  easier 
to  satisfy  as  shown  in  Fig.  5.20.  In  Fig.  5.20(a),  the  input  impedance  magnitude,  |Zin| , 

is  lower  than  the  magnitude  limit  only  at  f  =  103  Hz,  hence  the  phase,  ZZin ,  must  be 

checked  against  the  phase  specification  around  f  =  103  Hz.  As  shown  in  Fig.  5.20(b), 
the  phase,  ZZin ,  enters  slightly  the  unacceptable  phase  bands  only  at  point  A  where 

|Zin|  falls  below  the  magnitude  limit.  Therefore,  as  shown  in  Fig.  5.20(e),  the  polar 

plot  of  the  loop-gain  Tm  enters  slightly  the  forbidden  region  only  at  point  A, 


Fig.  5.19  Equivalent  DUT  input  impedance  specification  applied  to  the  conduction 
mode  of  Fig.  5.7(b),i.e.,  passive  DUT,  assuming  the  current  gain  is  function  of 
frequency,  i.e.,  k(s)  =  k0  (kjk12  / L2)/(s+  kjk12  /  L2),  where  ko=0.99.  Rp  is  2  Q. 
and  the  high  frequency  damping  circuit ,  Rp, ,  Cpl ,  has  not  been  added,  (a)  magnitude 
specification;  (b)  phase  specification;  (c)  Bode  magnitude  plot  of  Tm ;  (d)  Bode  phase 
plot  of  Tm ;  (e)  polar  plot  of  Tm  ;  and  (f)  zoom  of  (e)  around  the  origin  (on  the  next 
page). 
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Fig.  5.19  Continued,  (e)  polar  plot  of  Tm  ;  and  (f)  zoom  of  (e)  around  the  ori 


Fig.  5.20  Equivalent  DUT  input  impedance  specification  applied  to  the  conduction 
mode  of  Fig.  5.7(b),  i.e.  passive  DUT,  assuming  the  current  gain  is  a  function  of 
frequency,  i.e.,  k(s)  =  k0  (kjk^  /  L2)/(s  +  kjk,2  /  L2),  where  ko=0.99.  Rp  is  2  Q. 
and  the  high  frequency  damping  circuit  Rpl=lQ,  Cp,  =2uF  is  added,  (a) 
magnitude  specification;  (b)  phase  specification;  (c)  Bode  magnitude  plot  of  Tm ;  (d) 
Bode  phase  plot  of  Tm  ;  and  (e)  polar  plot  of  Tm  (on  the  next  page). 


Fig.  5.20  Continued,  (e)  polar  plot  of  T, 
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confirming  that  the  equivalent  DUT  input  impedance  almost  satisfies  the  input 
impedance  specification. 

For  the  active  dc-to-dc  regulator,  the  equivalent  small  signal  dynamic 
impedance  at  low-to-medium  frequency  will  be  negative  resistance  Rn  [60,61],  as 
shown  in  Fig.  5.7(c).  Figure  5.21  shows  the  equivalent  DUT  input  impedance 
specification,  without  the  addition  of  the  high  frequency  damping  circuit  Rpl,Cpl. 

This  input  impedance  specification  is  applied  to  the  negative  dynamic  load  resistance, 

Rn.  We  see  that  |Zin|  is  lower  than  the  magnitude  limit  around  frequencies  f  =  103 

Hz  and  f  =  lxl04  Hz  as  shown  in  Fig.  5.21(a).  The  phase  ZZin  must  then  be 
checked  against  the  phase  specification  around  these  frequencies.  As  shown  in  Fig. 
5.21(b),  the  phase  ZZin  enters  the  unacceptable  phase  bands  at  points  A  and  C  where 

|Zin|  falls  below  the  magnitude  limit.  The  polar  plot,  therefore,  of  the  loop-gain  Tm 
enters  the  forbidden  region  as  shown  in  Fig.  5.21(e),  although  the  integrated  system  is 
marginally  stable. 

If  the  high  frequency  damping  circuit  Rpi=lQ,Cpi=2uf  is  added  to  the  source 
subsystem,  however,  the  equivalent  DUT  input  impedance  specification  will  become 
easier  to  be  satisfied  as  shown  in  Fig.  5.22.  Figure  5.22(a)  shows  no  impedance 
overlapping,  therefore,  the  equivalent  DUT  input  impedance  will  satisfy  the  input 
impedance  specification.  As  shown  in  Fig.  5.22(e),  the  polar  plot  of  the  loop-gain  Tm 
always  remains  in  the  0.5  circle  and  does  not  enter  the  forbidden  region,  confirming 
that  the  equivalent  DUT  subsystem  satisfies  the  input  impedance  specification. 
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Fig.  5.21  Equivalent  DUT  input  impedance  specification  applied  to  the  conduction 
mode  of  Fig.  5.7(c),  i.e.,  negative  small  signal  dynamic  resistance,  assuming  the 
current  gain  is  a  function  of  frequency,  i.e.,  k(s)  =  k0  (kjk12  /  L2)/(s+  kjk12  /L2), 
where  k0  =  0.99 .  Rp  is  2  Q  and  the  high  frequency  damping  circuit,  Rpl ,  Cpl ,  has 

not  been  added,  (a)  magnitude  specification;  (b)  phase  specification;  (c)  Bode 
magnitude  plot  of  Tm ;  (d)  Bode  phase  plot  of  Tm  ;  (e)  polar  plot  of  Tm  ;  and  (f)  zoom 
of  (e)  around  the  origin  (on  the  next  page). 


Fig.  5.21  Continued,  (e)  polar  plot  of  Tm ;  and  (f)  zoom  of  (e)  around  the  origin. 
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Fig.  5.22  Equivalent  DUT  input  impedance  specification  applied  to  the  conduction 
mode  of  Fig.  5.7(c),i.e.  negative  small  signal  dynamic  resistance,  assuming  the  current 
gain  is  function  of  frequency,  i.e.,  k(s)  =  k0  (kjk12  /  L2)/(s+ kjk12  /  L2),  where 
k0  =  0.99 .  Rp  is  2  Q  and  the  high  frequency  damping  circuit,  Rpl  =  1Q,  Cpl  =  2uJF , 

is  added,  (a)  magnitude  specification;  (b)  phase  specification;  (c)  Bode  magnitude  plot 
of  Tm ;  (d)  Bode  phase  plot  of  Tm ;  and  (e)  polar  plot  of  Tm  (on  the  next  page). 


Fig.  5.22  Continued,  (e)  polar  plot  of  T, 
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5.6.4  Actual  PUT  Input  Impedance  Specification 

From  Fig.  5.15(b),  we  can  see  that  Zin  is  only  the  equivalent  DUT  impedance. 
The  actual  DUT  input  impedance  will  be 


Zjt  =  ZduL(1  _k) 
11  (1-k) 


where 


=  ZiB(l-k) 

=  ZDUT  (5.19) 
k  =  current  gain  of  the  controlled  current  source  (inverter) 


If  the  output  impedance  of  the  source  subsystem,  Zo,  is  multiplied  by  the  factor 
( 1  -k),  a  reduced  output  impedance  Zot  is  obtained 


Zot=Z0(l-k)  (5.20) 


The  loop-gain  Tm  of  the  integrated  system  can  be  defined  as 


(5.21) 


Therefore,  the  actual  DUT  input  impedance  specification  can  be  obtained  by 
substituting  Zot  for  Z0  in  Eq.  (5.17)  and  Eq.  (5.18). 
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If  the  equivalent  DUT  input  impedance  specification  in  Fig.  5.20,  which  has 
the  high  frequency  damping  circuit  added  and  is  applied  to  the  conduction  mode  of 
Fig.  5.7(b),  is  multiplied  by  (1-k),  the  actual  DUT  input  impedance  specification  can 
be  obtained  as  shown  in  Fig.  5.23  where  the  high  frequency  damping  circuit Rp,,Cp| 

is  also  added  and  the  conduction  mode  of  Fig.  5.7(b)  is  applied. 

Comparing  Fig.  5.20(a),  (b)  to  Fig.  5.23(a),  (b),  the  input  impedance 
specification  of  Fig.  5.23  is  only  an  impedance  level  shifting  of  that  in  Fig.  5.20  by  the 
factor  (1-k).  Since  k  *  ko  =  0.99  at  low-to-medium  frequency,  the  curves  in  these  two 

figures  will  be  similar  except  at  the  high  frequency  end,  i.e.,  for  f  >  104  Hz.  Through 
the  same  transformation,  using  formulas  from  Eq.  (5.19)  to  Eq.  (5.21),  the  equivalent 
DUT  input  impedance  specification  in  Fig.  5.22,  which  has  the  high  frequency 
damping  circuit  added  and  is  applied  to  the  conduction  mode  of  Fig.  5.7(c),  can  be 
transformed  into  the  actual  DUT  input  impedance  specification  as  shown  in  Fig.  5.24. 

- 

5.7  Discussion 

The  regenerative  load  system  is  transformed  into  two  subsystems,  the  utility 
source  subsystem  and  the  equivalent  DUT  subsystem.  The  stability  of  the  integrated 
system  can  be  analyzed  using  the  impedance  ratio  Tm=Zo/Zin  as  a  loop-gain  of  the 
system.    Hence,  the  impedance  overlap  may  cause  instability  of  the  system.  An 
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Fig.  5.23  DUT  input  impedance  specification  applied  to  the  conduction  mode  of  Fig. 
5.7(b),  i.e.  passive  DUT,  assuming  the  current  gain  is  function  of  frequency,  i.e., 
k(s)  =  k0(kik12 /L2)/(s+kik12 /L2),  where  ko=0.99.  Rp  is  2  Q  and  the  high 
frequency  damping  circuit,  Rpl=lQ,  Cpl=2uF,  is  added.  (a)  magnitude 
specification;  (b)  phase  specification;  (c)  Bode  magnitude  plot  of  Tm  ;  (d)  Bode  phase 
plot  of  Tm  ;  and  (e)  polar  plot  of  Tm  (on  the  next  page). 
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Fig.  5.23  Continued,  (e)  polar  plot  of  T, 
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Fig.  5.24  DUT  input  impedance  specification  applied  to  the  conduction  mode  of  Fig. 
5.7(c),i.e.  negative  small  signal  dynamic  resistance,  assuming  the  current  gain  is 
function  of  frequency,  i.e.,  k(s)  =  k0  (kjk12  /  L2)/(s+  kjk12  /  L2),  where  k0  =0.99. 
Rp  is  2  Q  and  the  high  frequency  damping  circuit  ,Rpl  =  1Q,  Cpl  =  2uF,  is  added, 
(a)  magnitude  specification;  (b)  phase  specification;  (c)  Bode  magnitude  plot  of  Tm ; 
(d)  Bode  phase  plot  of  Tm  ;  and  (e)  polar  plot  of  Tm  (on  the  next  page). 
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impedance  overlap  does  not,  however,  necessarily  imply  a  stability  problem.  In  this 
case,  the  Nyquist  criterion  can  be  applied  to  the  loop-gain  Tm  to  determine  the  system 
stability. 

If  the  load  were  passive,  the  phase  of  the  load  subsystem  would  always  be 

between  -90°  and  90° ,  which  is  generally  the  same  range  of  phase  that  the  output 
impedance  of  the  source  subsystem  has.  Therefore,  in  the  case  of  a  passive  load,  the 

phase  of  Tm  generally  cannot  be  worse  than  ±180° ,  so  the  polar  plot  of  the  loop-gain, 
Tm,  cannot  encircle  the  (-1,0)  point,  and  the  worst  possible  case  is  a  marginally 
stable  system.  However,  the  phase  characteristic  of  the  active  load  (e.g.  regulated  dc- 
to-dc  converter)  is  what  allows  the  interaction  between  the  source  and  load, 
subsystems  to  cause  potential  stability  problems  in  the  system.  For  the  active  load  as 
shown  in  Fig.  5.21,  Fig.  5.22,  and  Fig.  5.24,  the  phase  of  Tm  can  have  a  value  greater 

than  180°,  which  presents  the  opportunity  for  the  polar  plot  of  the  loop  gain  to 
encircle  the  (-1,0)  point  and  for  stability  problems  to  result. 

The  addition  of  the  damping  circuit  Rp,Cp  can  stabilize  the  system.  The  effect 
of  this  damping  circuit  is  to  stabilize  the  source  subsystem  and  at  the  same  time  to 
decrease  the  output  impedance  |Z0|  of  the  source  subsystem.  However,  due  to  the 
non-ideality  of  the  current  controller,  i.e.,  bandwidth  limitation,  a  second  impedance 

overlap  will  result  around  the  high  frequency  f  *  2  x  104  Hz  as  shown  in  Fig.  5.1 1(a), 
and  Fig.  5.13(a).  This  second  impedance  overlap  will  cause  the  system  to  become 
only  marginally  stable.  The  addition  of  a  high  frequency  damping  circuit  can  improve 
the  stability  remarkably,  as  shown  in  Fig.  5.16  and  Fig.  5.17. 
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After  the  high  frequency  damping  circuit  Rpi,Cpi  is  added  to  the  regenerative 
load  system,  the  magnitude  of  the  output  impedance  of  the  equivalent  source 
subsystem  will  be  decreased  to  almost  OdB,  i.e.,  1Q,  around  its  peak  value  in  the 

frequency  range  of  102Hz  to  5x  104Hz  as  seen  in  Fig.  5.16(a),  Fig.  5.17(a)  and  Fig. 
5.22(a).  For  the  actual  input  impedance  specification  of  the  DUT,  the  equivalent 
output  impedance  of  the  source  subsystem  is  at  the  level  of  -40dB ,  i.e.,  |Z0|  =  0.01 ,  as 
shown  in  Fig.  5.23(a)  and  Fig.  5.24(a).  This  means  that  for  the  stability  of  the  system, 
the  input  impedance  of  the  DUT  only  needs  to  be  larger  than  0.01Q,  i.e.,  |Zjt|  >  0.01, 

as  shown  in  Fig.  5.23(a)  and  Fig.  5.24(a). 

Additionally,  we  have  observed  that  in  Fig.  5.23(a)  and  Fig.  5.24(a)  an  input 

impedance  drop  of  the  DUT  around  frequency  f  =  10  Hz  exists.  This  is  due  to  the 
series  resonance  of  the  leakage  inductance  of  the  transformer  and  the  large  capacitance 
of  C  in  Fig.  5.7(b)  and  (c).  Thus,  the  input  impedance  of  DUT  will  be  approximately 
equal  to  series  resistance  of  the  transformer  and  equivalent  core  loss  rft  where 
r/j  =0.001  and  r/t  =0.01  as  assumed  in  the  analysis.  However,  in  a  real  situation,  the 
coil  resistance  of  the  transformer  will  be  greater  than  0.00 1Q  and,  together  with  the 
equivalent  resistance  of  the  core  loss  r/t,  will  make  the  |Zit|  larger  than  0.1Q  at  this 

resonance  frequency  f*103Hz.  Therefore,  the  stability  of  this  regenerative  load 
system  will  be  guaranteed  for  all  practical  DUT. 


5.8  Simulation  Results 
Figure  5.25  shows  the  PSpice  simulation  results,  when  the  low-to-medium 
frequency  damping  circuit,  Rp=10Q  Cp=10uF,  and  the  high  frequency  damping 
circuit,  Rpi  =  1Q  Cp,  =  2uF ,  are  added  to  the  regenerative  load  system  with  source 

inductance  Ls  =  2mH.  From  Fig.  5.25,  we  saw  that  all  the  waveforms  are  clean,  i.e., 
there  is  no  potential  instability  in  the  system.  Finally,  the  transient  response  is  fast  and 
good,  i.e.,  less  than  two  cycles  after  the  steady-state  of  DUT. 

5.9  Experimental  Results 
Figure  5.26  shows  the  simulation  results,  where  Vs  =  32  v  RMS  (voltage 

waveform  pre-scaled  down  by  1/3),  the  source  inductance  Ls=  the  leakage  inductance 
of  power  transformer,  the  low-to-medium  frequency  damping  circuit  Rp=3.3Q 
Cp=12uF,  and  the  high  frequency  damping  circuit  Rpl=2Q   Cpl=2|iF  and 

Cdc=680uF.  The  energy-storage  capacitor  voltage  is  set  at  100V.  The  DUT  is  as  that 

in  Fig.  4.10.  The  output  voltage  of  the  DUT  is  40  V  and  the  output  current  of  the  DUT 
began  with  a  value  of  0.3A,  rose  up  to  2A  and  fell  back  again  to  0.3A.  From  Fig.  5.26, 
we  see  that  all  the  waveforms  are  clean,  i.e.,  there  is  no  potential  instability  in  the 
system. 
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Fig.  5.25  Simulated  current  and  voltage  waveforms  for  the  regenerative  load  system  with 
utility  source  inductance  Ls  =2mh  ,  after  both  the  low-to-medium  frequency  damping 
circuit  "Rp=10Q  Cp=10uT"  and  high  frequency  damping  circuit  "Rpl=lQ  Cpl=2uF" 
are  added. 

(a)  Utility  source  voltage  vs  =  v(7,3) ;  utility  source  current  is  =  i(vms) 

(b)  DUT  input  current  i  L  =  i(vm) 

(c)  Current  controller  current  icc  =  i(vmc) 

(d)  DUT  output  voltage  v0  =  v(l) ;  DUT  output  current  i0  =i(vmle) 

(e)  Capacitor  (Cdc)  voltage  Vdc  =  v(71,FGnd) 

DUT  output  current  changes  from  a  peak  value  of  OA  to  10A  at  the  fifth  cycle  and  from 
1  OA  to  OA  at  the  tenth  cycle. 
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Fig.  5.26  Experimental  results  for  the  regenerative  load  system  applied  to  the  conduction 
mode  of  Fig.  5.7(b),  where  Cc=  680uF,  Ls=  transformer  leakage  inductance,  Cp=12uF, 
Rp=3.3Q,  Cp,=2uF,  Rp,=2Q; 

(a)  trace  1;  utility  source  voltage  vs  =  v(7,3),  (pre-scaled  down  by  1/3) 
trace  2; ;  utility  source  current  is  =  i(vms),  (current  scale  =  1(A/V)) 
trace  3:  DUT  input  current  iL  =  i(vm),  (current  scale  =  1(A/V)) 

trace  4:  current  controller  current  icc  =  i(vmc).  (current  scale  =  1(A/V)) 

(b)  trace  1:  utility  source  voltage  vs  -  v(7.3),  (pre-scaled  down  by  1/3) 
trace  2:  DUT  output  voltage  v0  =  v(l) 

trace  3:  DUT  output  current  i0  =  i(vmle),  (current  scale  =  l(A/V)) 
trace  4;  energy-storage-capacitor  (Cdc)  voltage  Vdc  =  v(71,FGnd) 


147 

5.10  Summary 

In  this  chapter,  we  presented  the  input  impedance  specification  of  DUT  for  the 
system  stability.  At  first,  the  whole  regenerative  system  is  transformed  into  two 
subsystems,  using  the  proposed  equivalent  impedance  concept.  This  equivalent 
impedance  concept  combines  the  controlled  current  source  with  its  associated 

Z 

impedance  to  become  an  equivalent  impedance,  Zea  =  ,  where  Z  is  the  original 

impedance  and  k  is  the  gain  of  controlled  current  source.    After  this  equivalent 
impedance  transformation,  the  whole  system  can  be  divided  into  two  equivalent 
subsystems.  These  two  subsystems  are  equivalent  source  subsystem  and  equivalent 
DUT  subsystem. 

The  two  subsystem  interaction  theorem  [57,66]  is  employed  to  analyze  the 
stability  of  the  integrated  system  using  the  impedance  ratio  as  the  equivalent  loop-gain 
of  the  system.  It  is  found  that  for  the  passive  DUT  load,  the  phase  of  the  equivalent 
2 

loop-gain,  Tm  =  — — ,  cannot  be  worse  than  ±180° ,  so  the  polar  plot  of  the  loop-gain 
^in 

Tm  cannot  encircle  the  (-1,0)  point,  and  the  worst  possible  case  is  a  marginally  stable 
regenerative  load  system.  However,  the  phase  characteristic  of  the  small  signal 
dynamic  input  impedance  of  the  active  DUT  (e.g.  regulated  dc-to-dc  converters)  is 
what  allows  the  interaction  between  the  source  and  DUT  subsystems  to  cause  potential 
stability  problems  in  the  system.  Additionally,  the  bandwidth  limitation  of  the  current 
controller  has  an  adverse  effect  on  the  system  stability. 
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We  found  that  the  addition  of  a  second  high  frequency  damping  circuit  to  the 
system  can  increase  the  stability  margin  of  the  system.  This  two-damping-circuits 
concept  can  also  be  applied  to  other  places  of  circuit  design  and  is  one  of  the 
contributions  of  this  study. 

With  the  addition  of  these  two  damping  circuits,  low-to-medium  frequency 
damping  circuit  and  high  frequency  damping  circuit,  the  proposed  regenerative  load 
system  will  have  enough  stability  margin  for  the  practical  DUT,  even  for  the  case  of 
active  power  regulators.  The  simple  and  enough  condition  for  the  system  stability  is 
that  the  magnitude  of  the  input  impedance  of  the  DUT  needs  to  be  greater  than  0.1Q 
for  all  the  frequency  range  of  interest. 


CHAPTER  6 


ACTIVE  DAMPING  AND 
COMPENSATION  OF 
THE  REACTIVE  CURRENT 

6.1  Introduction 

In  Chapter  5,  the  two  damping  circuits  were  needed  to  make  the  regenerative 
load  system  satisfy  the  DUT  input  impedance  specification.  These  two  damping 
circuits  are  Rp,Cp  and  Rpi,Cpi  Since  these  circuits  are  passive,  they  will  cause 
energy  dissipation  in  the  resistors  Rp  and  Rpi  and  decrease  the  total  efficiency  of  the 
system.  Also,  these  damping  circuits  can  absorb  the  capacitive  reactive  current  and 
cause  the  power  factor  of  the  source  current  is  to  be  less  than  one,  i.e.,  P.F.  <  1 .0. 

The  active  damping  circuit  is  introduced  in  Sec.  6.2  to  replace  the  passive 
damping  circuits.  The  energy  consumption  in  the  active  damping  resistor  Rp  act  is 

restored  in  the  energy  storage  capacitor  Cdc .  Therefore,  the  efficiency  of  the  system 
can  be  increased. 

For  the  compensation  of  the  reactive  current  in  the  damping  circuits,  the 
current  controller  needs  to  supply  additional  compensation  current  icomp-60Hz  t0 

balance  out  the  reactive  component  current  produced  by  the  damping  circuits  at  the 
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frequency  of  60  Hz.  The  power  factor  can  then  be  increased  to  nearly  one,  i.e.,  P.F.  * 
1.0.  The  details  are  presented  in  Sec.  6.3. 

The  switching  frequency  of  the  current  controller  is  fs  =  20kHz  in  this  design. 

The  bandwidth  of  the  current  controller  will,  therefore,  be  smaller  than  -  fs  according 

to  the  sampling  theorem  [72].  In  practice,  however,  the  bandwidth  is  further  limited  to 

approximately  —  fs  [69,70].   This  bandwidth  limitation  will  restrict  the  frequency 

range  of  the  active  damping  circuit  to  the  low-to-medium  frequency.  Hence,  only  the 
low-to-medium  frequency  damping  circuit  Rp,Cp  can  be  replaced  with  the  active 
damping  circuit.  For  the  high  frequency  damping  circuit  Rpi.Cpi  the  passive  damping 
is  still  needed.  The  bandwidth  of  the  current  controller  is  discussed  Sec.  6.4.  The 
whole  composite  control  circuit  for  the  active  damping  and  compensation  of  the 
reactive  current  is  presented  in  Sec.  6.5. 

The  validity  of  this  modified  current  controller  design  with  active  damping  and 
reactive  current  compensation  is  verified  by  using  a  PSpice  simulation.  The  results 
show  that  the  regenerative  load  system  is  still  stable  and  has  better  overall  efficiency 
and  nearly  unity  power  factor,  i.e.,  P.F.  «  1.0. 

6.2  Active  Damping 
Since  the  resistor  Rp  in  the  passive  damping  circuit  Rp,Cp  in  Fig.  6.1(a)  will 
dissipate  energy,  the  passive  damping  circuit  can  be  replaced  with  an  active  damping 
circuit  as  shown  in  Fig.  6.1(a),  represented  by  the  current  source  Ipacl. 
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Fig.  6. 1  The  active  damping  circuit. 
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The  current  in  the  active  damping  circuit,  Ipact  must  be  equal  to  Ippassive  of  the 
passive  damping  circuit. 


Ip.act     1  p.  passive 


(6.1) 


Therefore,  the  current  controller  needs  to  provide  additional  current,  Ip-act.  The 


reference  control  current  I*ef  for  the  current  controller  will  be 


^ref  —  ^i,scale^(^cl  +^in)  ^ref.damp 


(6.2) 


where  the  minus  sign  comes  from  the  direction  of  Ipact  opposite  to  the  direction  of 
k(Ici  +Ijn)-  The  additional  reference  control  current  Iref  damp  of  the  active  damping 


will  be 


T  —  V  T 

1ref,damp     "-i.scale1  p.active 


-  k  T 

"N.scale1?,  passive 


(6.3) 


where 


kj  scale  =  current  scale  of  the  current  controller 


Hence,  the  I  p,passive  can  be  calculated  as  follows 
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where 


'  p, passive 


^v, scale    .  , 
7  k  v, scale ^z.scale 

^pKz,scale 


v(81,0) 


~    7         ^v, scale  ^z.scale  (6.4) 
^p,eq 


Zp,eq  —  ^p.eq  +   p  ^z,scale^p 
SUp,eq 


v 

kv       =  5 —  =  voltage  scale  in  Fig.  6.1(b) 

v,scale  v(g1>0) 


:z,scaie  =         =  impedance  scale  in  Fig.  6.1(b) 
Zp 


Since  v^'^  =  v^^'^  as  in  Fig.  6.1(b),  Eq.  (6.4)  can  also  be  written  as 
Zp,eq  Rpeq 


'p, passive      p         *v,  scale  *-z,scale 
p.eq 

v(82,0)  _ 

K  v, scale 


lvp,eq  '  ^-z.scale 


v(82;0) .  6 

D       *-v,scale  IP'?) 


Substituting  Eq.  (6.5)  into  Eq.  (6.3)  results  in 


^ref.damp  —  v(82,0)k  v  scaje  ^i.scale 


P 


where 


=  v(82,0)Av  (6.6) 


A   -k      ,    1  1 

v  v.scale 


R  k 

scale 

For  this  design,  the  following  scaling  are  assumed 

v         — 10 

v.scale  ~  JU 

ki.sca.e=0.1  [V/A] 

and  the  damping  resistor  Rp  is  equal  to  10Q.  Hence,  Av  in  Eq.  (6.6)  will  be  0.3.  In 
addition,  the  relation  between  Zp  and  Zp,eq  leads  to 

Rp.eq  =  kz,scaleRp  (6.7) 
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6.3  Compensation  of  the  Reactive  Current 
In  order  to  balance  out  the  reactive  current  component  of  the  damping  circuit, 
Zp  and  Zpi,  the  current  controller  must  supply  another  60  Hz  compensation  current, 

1comp,60Hz ' 

such  that  icomp,60Hz  =  idamp.60Hz .  as  shown  in  Fig-  6-4-  Then  the  Power 

factor  of  the  utility  source  current,  is,  will  approach  one,  i.e.,  P.F.w  1.0. 

Figure  6.2  shows  the  calculation  circuit  for  the  reference  60  Hz  compensation 

current,  I  refcomp,  which  is  added  to  the  reference  control  current  I*ef  of  the  current 

controller.  The  main  concern  of  Fig.  6.2  is  that  the  voltage  across  the  equivalent 
impedance  Zpeq  and  Zpi  eq  must  be  60  Hz  pure  sinusoidal  waveform  and  in  phase 

with  the  utility  voltage,  vs.  This  voltage,  v(59a2,0)  in  Fig.  6.2(a),  is  obtained  by 
multiplying  the  phase  locked  signal  vsr  with  the  scaled  down  peak  value  of  the  source 
voltage,  v(59al,0),  as  shown  in  Fig.  6.2(c).  The  low  pass  filter  in  Fig.  6.2(c)  is  used  to 
filter  out  the  high  frequency  oscillation  noise.  Therefore,  the  voltage  v(59a2,0)  will  be 
a  60  Hz  pure  sinusoidal  voltage  and  in  phase  with  the  voltage  source,  vs .  The  voltage 

scale  between  vs  and  v(59a2,0)  is  kv  scale  =  vs/v(59a2,0)  =  30  for  this  design. 

Through  a  similar  derivation  process  as  in  Sec.  6.2,  the  reference  60  Hz 
compensation  current  for  the  current  controller  is  obtained  as 
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Fig.  6.2  The  calculation  circuit  for  the  reference  compensation  current  I 

(a)  Block  diagram 

(b)  Sampling  control  logic 

(c)  Full-wave  rectifier  and  low  pass  filter 
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Fig.  6.2  The  calculation  circuit  for  the  reference  compensation  current  ^  . 

(b)  Sampling  control  logic 

(c)  Full-wave  rectifier  and  low  pass  filter 
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1  ref ,zp.60Hz  -  v(85,0)k  v  sca)e  — —  k ,  sca[e  (6.9) 

p 

=  v(86,0)  in  Fig.  6.2(a) 


\  ref  ,zpl.60Hz  =  v(87,0)k  v  scale  — —  k  j  scaie  (6. 1 0) 

Kpl 

=  v(88,0)  in  Fig.  6.2(a) 

where 

V 

kv  sca)e  =  voltage  scale  =  5  =30  used  in  this  design 

v(59a2,0) 

kj  scale  =  current  scale  =  ref'comp  =  0.1(v/a)  used  in  this  design 

^comp 


Also,  the  relation  between  Zp  and  Zpl  eq  is  as  follows 

^pl.eq  =  ^z,scaleRpl  (6.1 1) 

1 


^pl,eq^7  Cpl  (612) 

K  z.scale 

scale  =    P'  eq  =  impedance  scale 
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6.4  Bandwidth  of  the  Current  Controller 
From  Eq.  (4.6),  it  can  be  seen  that  the  bandwidth  of  the  current  controller  could 
be  as  high  as 


fB  =  ^-^r^2-  =  30kHz  (6.13) 


2n  L 


2 


However,  according  to  the  sampling  theorem  [72],  the  bandwidth  of  the  current 
controller  will  be  limited  to  half  of  its  switching  frequency  as 

fB<^fs  =  10kHz  (6.14) 

where  fs  =  20kHz  is  the  PWM  switching  frequency  . 

In  practice,  the  bandwidth  of  the  current  controller  is  further  reduced  in  Eq. 
(6.15)  [69] 

fB*Ifs  (6.15) 


Figure  6.3  shows  the  PSpice  simulation  results  of  the  current  controller.  In 
Fig.  6.3(a),  the  frequency  of  the  control  current  I*ef  is  set  at  f  =  2kHz,  where  we  noted 
that  the  output  current  waveform  of  the  current  controller  becomes  triangular  and 
distorted.  If  the  frequency  is  decreased  to  f  =  1 .5  kHz,  the  output  current  waveform  in 
Fig.  6.3(b)  still  shows  notable  distortion. 


Fig.  6.3  Simulation  results  of  the  bandwidth  of  the  current  controller. 

(a)  Reference  control  current  I*ef  at  f  =  2  kHz 

(b)  Output  current  of  current  controller  at  f  =  2  kHz 

(c)  Reference  control  current  I*ef  at  f  =  1.5  kHz 

(d)  Output  current  of  current  controller  at  f  =  1 .5  kHz 

(e)  Reference  control  current  I*ef  at  f  =  1  kHz 

(f)  Output  current  of  current  controller  at  f  =  1  kHz 
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When  the  frequency  is  further  reduced  to  f  =  1kHz,  the  result  in  Fig.  6.3(c) 
shows  that  the  output  current  waveform  is  almost  sinusoidal  and  the  distortion  is 
small.  Therefore,  the  actual  bandwidth  of  this  current  controller  will  be  slightly  higher 
than  1  kHz,  based  on  the  distortion  criterion 

fB*UkHz  (6.16) 

This  bandwidth  limitation  of  the  current  controller  will  restrict  the  emulation  of 
the  active  damping  circuit  to  only  the  low-to-medium  frequency  damping,  i.e.,  Rp,Cp 
circuit.  However,  for  the  high  frequency  damping,  the  passive  damping  circuit 
Rpi  Cpi  is  still  needed,  as  described  in  the  following  section. 

6.5  Combination  of  the  Active  Damping  and 
Compensation  of  the  Reactive  Current 

The  active  damping  circuit  and  the  compensation  of  reactive  current  can  be 
combined  together  and  built  into  the  current  controller  as  shown  in  Fig.  6.4.  Because 
of  the  bandwidth  limitation  of  the  current  controller,  fB  W  «  1.1  kHz,  only  the  low-to- 
medium  frequency  damping  circuit  Rp,Cp  can  be  built  as  an  active  damping  circuit. 
The  high  frequency  damping  circuit  Rp^Cpi  still  needs  to  remain  as  a  passive  circuit. 
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Fig.  6.4  The  incorporation  of  active  damping  and  compensation  of  reactive 
current  into  the  current  controller. 

(a)  main  power  circuit 

(b)  calcuation  circuit  for  the  reference  control  current  I,.  * 

ret 
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Figure  6.4(a)  shows  the  block  diagram  of  the  main  power  circuit  of  the 
regenerative  load  system.  The  calculation  circuit  for  the  reference  control  current  I*ef , 
is  shown  in  Fig.  6.4(b)  where  I*ef  consists  of  three  components. 


^ref  —  ki,scalekOcl  +  Ijn)  +  ^  ref  ,comp,60Hz     1  ref  ,damp,zp  (6-17) 


where 


ki.scale  =  current  scale  of  the  current  controller 
The  minus  sign  in  front  of  Iref^amp.zp  's  due  t0  tne  reference  direction  of  the 
Iref  damp  zp  and  is  downward,  opposite  to  the  reference  direction  of  the  controller. 

The  value  of  the  equivalent  impedance  Zpeq  in  Sec.  6.3  is  determined  as 

follows 


Zp,eq  ~  Rp.eq  +  ~        ~  kz,sca|e (Rp  +  ~    )  (6.18) 
S^p.eq  SCp 


For  kz,scaie  =  100>  men  RP.eq  and  Cpeq  are 


Rpeq  =  100Rp  =  1  kQ 


C^-^Cp-O-liiF  (6.19) 
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Similarly,  using  kzsca,e  =  10,  the  Rpleq  and  CpUq  in  the  compensation  current 
calculation  circuit  are 

Rpleq=10Rp=  ion 

Cp,,eq=^Cp  =  0.2uF  (6.20) 

However,  due  to  the  bandwidth  limitation  of  the  current  controller,  the  actual  value  of 
cP,eq  adopted  is 

Cp,eq  =  0.2uF             Cp  =  20uF  (6.21) 
Therefore,  the  3dB  frequency  of  the  active  damping  circuit  will  be 


f  _ 


3dB    27iRpeqCpeq  2nRpCp 

-  796  Hz  (6.22) 

This  f3dB  frequency  is  within  the  bandwidth  of  the  current  controller  . 

This  composite  circuit  has  two  advantages:  (1)  a  nearly  unit  power  factor,  P.F.« 
1.0,  of  the  source  current,  is,  due  to  the  compensation  of  the  reactive  current  and  (2) 
low  energy  dissipation  in  the  damping  circuit,  because  the  only  energy  dissipation  in 
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the  passive  damping  circuit  Rpi  Cpi  is  less  than  0.1  W.  These  advantages  are  verified 
through  PSpice  simulation. 

6.6  Simulation  Results 
Figure  6.5  shows  the  simulation  results,  using  PSpice  software,  where  the 
active  damping  circuit  Cp=20uF,  Rp=10Q,  high  frequency  passive  damping  circuit 
Cpi=2uF,  Rpi=lQ,  with  the  reactive  current  compensation.  The  value  of  capacitance 
Cp  is  changed  from  the  original  lOuF  to  20uF  because  of  the  bandwidth  limitation  of 
the  current  controller.  If  the  Cp=10uF  is  used,  the  f3dB  of  the  active  damping  circuit 
will  be 

f  1       -  1  

3dB     27TRpCp  27r.xl0xl0xl(T6 

s  1.6  kHz  (6.23) 

This  frequency  is  higher  than  the  actual  bandwidth  of  the  current  controller,  fBw  *  1.1 
kHz  as  in  Eq.  (6.16).  Hence,  the  PSpice  simulation  shows  an  instability  phenomenon. 
If  Cp=20uF  is  used,  the  f3dB  will  be  approximately  0.8  kHz,  within  in  the  bandwidth 
of  the  current  controller.     However,  the  replacement  of  Cp=10uF  with  20uF 
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Fig.  6.5  Simulated  results  of  the  regenerative  load  system  with  the  active  damping  and 
compensation  of  reactive  current  components.  The  DUT  is  non-regulated. 

(a)  Utility  source  voltage  vs  =  v(7,3) ;  utility  source  current  is  =  i(vms) 

(b)  DUT  input  current  iL  =  i(vm) 

(c)  Current  controller  current  icc  =  i(vmc) 

(d)  DUT  output  voltage  v0  =  v(l) ;  DUT  output  current  i0  =  i(vmle) 

(e)  Capacitor  (Cdc)  voltage  Vdc  =  v(71,FGnd) 

DUT  output  current  changes  from  a  peak  value  of  OA  to  10A  at  the  fifth  cycle  and  from 
10A  to  OA  at  the  tenth  cycle. 
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is  equivalent  to  the  replacement  of  Cp  eq  =  0.1  uF  with  0.2uF.  The  cost  of  the  circuit 
almost  remains  unchanged.  The  simulation  results  of  Fig.  6.5  show  that  the  reactive 
current  component  of  utility  source  current,  is,  is  very  small  and  the  power  factor  is 
nearly  one,  i.e.,  P.F.«  1.0.  In  addition,  the  power  dissipation  in  the  passive  high 
frequency  damping  circuit  is  less  than  0. 1  W. 

In  Fig.  6.6,  the  DUT  is  replaced  with  a  regulated  Buck  converter.  This  is  seen 
from  the  output  voltage  v(5)  in  the  simulation  of  the  Buck  converter,  which  remains 
constant  during  the  transient  response.  Also,  the  regenerative  load  system  is  stable, 
even  though  the  dynamic  small  signal  resistance  of  the  regulated  Buck  converter  is 
negative  in  the  low  -to-medium  frequency  range  [60,61]. 

6.7  Summary 

Because  the  damping  circuit  is  used  to  stabilize  the  regenerative  load  system, 
these  conventional  passive  damping  circuits  will  have  power  dissipation  and  also 
produce  line  frequency  (60  Hz)  reactive  current.  Therefore,  the  overall  efficiency  of 
the  system  will  be  decreased  and  the  power  factor  of  the  utility  source  current,  is,  is 
reduced.  In  this  chapter,  an  active  damping  technique  was  proposed,  using  the 
impedance  emulator  concept  to  absorb  the  power  into  the  energy  storage  capacitor  Cdc 
and  reuse  it.  Hence,  the  overall  efficiency  of  the  system  can  be  increased.  Owing  to 
the  bandwidth  limitation  of  the  current  controller,  only  the  low-to-medium  frequency 
damping  circuit  can  be  replaced  with  the  active  damping  circuit.    However,  the 
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Fig.  6.6  Simulated  results  of  the  regenerative  load  system  with  the  active  damping  and 
compensation  of  reactive  current  components.  However  the  DUT  is  regulated  "Buck" 
converter. 

(a)  Utility  source  voltage  vs  =  v(7,3) ;  utility  source  current  is  =  i(vms) 

(b)  DUT  input  current  iL  =  i(vm) 

(c)  Current  controller  current  icc  =  i(vmc) 

(d)  DUT  output  voltage  v0  =  v(  1 );  DUT  output  current  i0  =i(vmle) 

(e)  Capacitor  (Cdc)  voltage  Vdc  =  v(71,FGnd) 

DUT  output  current  changes  from  a  peak  value  of  OA  to  10A  at  the  fifth  cycle  and  from 
10A  to  OA  at  the  tenth  cycle. 


remaining  passive  high  frequency  damping  circuit  has  power  dissipation  less  than  0.1 
W. 

In  addition,  the  line  frequency  (60  Hz)  reactive  current  in  the  damping  circuit 
can  also  be  compensated  by  the  compensation  current  supplied  from  the  current 
controller.  Both  the  active  damping  and  the  compensation  of  the  line  frequency 
reactive  current  can  be  built  into  the  current  controller,  using  the  superposition 
principle.  The  merits  of  this  composite  design  are  (1)  an  increase  in  the  overall  system 
efficiency  and  (2)  an  increase  of  the  power  factor  of  the  utility  source  current  to  nearly 
one.  The  feasibility  of  this  scheme  is  verified  by  PSpice  simulation.  The  results 
demonstrate  that  the  integrated  system  is  stable  even  during  the  transient  response  and 
nearly  unity  power  factor  of  the  utility  source  current  is  obtained.  Also,  the  overall 
efficiency  is  improved. 


CHAPTER  7 
CONCLUSIONS  AND  FUTURE  WORKS 

7.1  Conclusions 

This  dissertation  has  presented  the  analysis  and  design  of  the  regenerative  load 
system.  This  system  consists  of  two  parts:  (1)  an  impedance  emulator  and  (2)  an 
inverter  with  power  filtering  function.  The  proposed  system  can  be  used  for  long-term 
testing  of  power  equipment  and  motor-generator  set  and  has  the  following  merits. 

1.  Economic,  because  of  the  recycle  of  electrical  energy,  the  energy  consumption  will 
be  less  than  20%  of  the  conventional  load.  If  the  soft  switching  scheme  is 
employed  for  the  converters,  the  energy  dissipation  can  be  even  lower  than  10%. 

2.  Safety,  due  to  the  low  energy  dissipation,  no  appreciable  Joule  heating  is  observed. 

3.  Power  quality,  since  the  current  of  the  utility  source  is  shaped  into  pure  sinusoidal 
waveform  with  p.f.  «  1.0,  problems  such  as  harmonic  currents  and  EMI  will  be 
greatly  reduced. 

In  Chapter  2,  a  new  APF  technique  is  proposed,  using  integration  of  power  and 
a  sampling  technique  to  simplify  the  calculation  of  the  real  fundamental  component  of 
the  load  current.  Also,  the  energy-balance  concept  in  the  energy-storage  capacitor  is 
used  to  simplify  the  design  of  the  conventional  APF-capacitor-voltage-control  circuit. 
The  merits  of  this  APF  are  (1)  the  calculation  of  the  current  is(t)  required  by  the 
utility  source  is  simplified;  (2)  a  larger  ripple  voltage  can  be  tolerated  in  the  energy 
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storage  capacitor,  therefore,  a  smaller  energy  storage  capacitor  is  needed;  and  (3)  only 
a  proportional  capacitor-voltage-control  method  is  applied  in  the  controller  design. 
Hence,  the  transient  response  is  fast  and  good  (less  than  two  cycles  of  60  Hz). 

For  the  stability  analysis  of  the  regenerative  load  system,  we  found  that  with 
the  existence  of  utility  source  inductance,  Ls,  the  system  will  become  unstable.  After 
the  addition  of  a  damping  circuit,  the  system  can  be  stabilized.  However,  even  with 
the  addition  of  a  damping  circuit,  the  current  gain  of  the  APF  must  be  kept  below  one, 
i.e.,  k  <  1.0,  for  the  stability  of  the  regenerative  load  system.  The  reason  for  this 
phenomenon  is  that  if  the  current  gain  of  the  current  controller  is  larger  than  one,  the 
system  will  be  a  positive-feedback  system  with  loop-gain  greater  than  one.  Under  this 
condition,  this  system  is  unstable.  This  conclusion  has  been  verified  by  using  the  root- 
locus  analysis  and  PSpice  simulation.  Therefore,  the  criterion  for  a  stable  APF  or 
regenerative  load  system  is  that  the  damping  circuit  must  be  added,  and  at  the  same 
time,  the  current  gain  of  the  current  controller  needs  to  be  smaller  than  one,  i.e.,  k  <  1 . 

In  order  to  generate  the  input  impedance  specification  of  the  DUT  for  the 
system  stability,  the  whole  system  is  transformed  into  two  subsystems  using  the 
proposed  equivalent  circuit  concept.  This  equivalent  circuit  concept  is  to  merge  the 
controlled  current  source  into  the  associated  impedance  to  form  an  equivalent 
2 

impedance,  Zeq  =  ,  where  Z  is  the  original  impedance  and  k  is  the  gain  of  the 

(1  —  k) 

controlled  current  source.  Therefore,  the  system  can  be  divided  easily  into  two 
subsystems.  Thus  the  two  subsystem  interaction  theorem  [57,66]  can  be  employed  to 
analyze  the  system  stability.  In  Chapter  5,  a  methodology  for  deriving  the  DUT  input 
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impedance  specification  is  presented.  We  also  found  that  the  addition  of  a  second  high 
frequency  damping  circuit  can  increase  the  stability  margin.  This  two-damping- 
circuits  concept  is  a  new  contribution  of  this  work  and  can  be  applied  to  other  places 
of  circuit  design.  With  these  two  damping  circuits  added,  it  is  guaranteed  that  the 
whole  system  will  be  stable  for  all  practical  DUT. 

The  conventional  passive  damping  circuits  will  consume  energy  and  decrease 
the  overall  system  efficiency.  In  addition,  the  passive  damping  circuits  can  absorb  the 
reactive  current  and  degrade  the  power  factor  of  utility  current.  Therefore,  in  Chapter 
6,  an  active  damping  circuit  is  proposed  to  replace  the  passive  one.  Due  to  the 
bandwidth  limitation  of  the  controlled  current  source,  only  the  low-to-medium 
frequency  damping  circuit  can  be  replaced  with  the  active  damping  circuit.  The  high 
frequency  damping  circuit  is  still  passive.  However,  the  power  dissipation  of  this 
damping  circuit  is  less  than  0.1W.  In  addition,  the  reactive  current  of  the  damping 
circuits  is  compensated  by  the  compensation  current  supplied  from  the  current 
controller.  The  merits  of  this  composite  current  controller  design  are  (1)  the  increase 
of  the  system  efficiency  and  (2)  the  power  factor  of  the  utility  source  can  be  increased 
to  nearly  one,  i.e.,  p.f.  *  1.0. 

Further  research  is  recommended  in  the  following  areas: 

•  Design  of  the  bandwidth  of  the  Active  Power  Filter  and  analysis  of  the  system 
stability  using  more  detailed  current-controller  model. 

•  Investigation  of  the  pole  placement  as  an  alternative  analysis  and  design  method. 
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